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We derive general closed-form analytical expressions for the finite-energy one- and two-electron 
spectral-weight distributions of an one-dimensional correlated metal with on-site electronic repulsion. 
Our results also provide general expressions for the one- and two-atom spectral functions of a 
correlated quantum system of cold fermionic atoms in a one-dimensional optical lattice with on- 
site atomic repulsion. In the limit of zero spin density our spectral-function expressions provide 
the correct zero-spin density results. Our results reveal the dominant non-perturbative microscopic 
many-particle mechanisms behind the exotic spectral properties observed in quasi-one-dimensional 
metals and correlated systems of cold fermionic atoms in one- dimensional optical lattices. 

PACS numbers: 70 

I. INTRODUCTION 

The main goal of this paper is to provide a detailed derivation of the one-electron and two-electron spectral-weight 
distributions recently used in Refs. [lj, |2|, |3j in the unusual spectral properties of low-dimensional metals. Moreover, 
our results will be used elsewhere in the study of the spectral properties of the new quantum systems described by 
cold fermionic atoms in optical lattices: following the experimental studies of strongly correlated quantum systems 
of ultra cold bosonic atoms Q, new experiments involving cold fermionic atoms (such as 6 Li) in an optical lattice 
formed by interfering laser fields are in progress |]|. 

We profit from symmetries and properties specific to the metallic phase of the ID Hubbard model |fj to derive 
closed-form analytical expressions for the finite-energy one- and two-electron spectral- weight distributions. The model 
also describes cold fermionic atoms in optical lattices, provided that the electrons are replaced by these atoms. The 
recently developed quantum optical tools for manipulating atoms can be used to realize a ID lattice or any other lattice 
described by the Hubbard model 0, 0, 13 ■ Such correlated quantum systems have applications in quantum computing 
0- Our study involves the use of the holon and spinon description of Ref. and of the general pseudofermion 
dynamical theory introduced in Ref. [lCj : we are able to explicitly calculate the pseudfermion spectral weights for the 
metallic phase. The pseudofermion description was introduced in Ref. 0] an d is such that the pseudofermions which 
describe the elementary excitations are freely propagating, and the information about the interactions is encoded in 
their overall phase shifts The method used here is a generalization for all values of the on-site repulsion U of 
that used in Refs. 0,0] for U — ► oo. For finite values of U there were no closed-form analytical expressions for the 
finite-energy spectral- weight distributions, in contrast to simpler models 0. 

Our general closed-form analytical expressions for the one- and two-electron spectral functions are controlled by the 
dominant microscopic processes which are behind nearly the whole spectral weight for all values of the momentum 
k, energy u>, and on-site repulsion U. In addition to other spectral features, we explicitly derive the spectral-weight 
distribution expressions in the vicinity of the the singular and edge branch lines used in the study of the photoemission 
spectra and phase diagram of quasi-lD metals in Refs. [HI3.l3.1l^. The preliminary results of those papers confirm 
that the spectral-weight distributions derived here describe quantitatively, for the whole energy band width, the 
unusual spectral properties observed in such materials. 

The paper is organized as follows: In Sec. [HI we introduce the model and provide basic information about the 
pseudofermion description needed for our studies. The specific classes of elementary microscopic processes which 
generate the one and two-electron spectral weight of the metallic phase is the subject of Sec. IIIII In Sec. IIVI 
we calculate the general expressions for the finite-energy spectral-weight distributions of the metallic phase. The 
characterization of the dominant processes behind such distributions and the derivation of the closed-form analytical 
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expressions generated by these processes for the same distributions is the subject of Sec[V] Finally, the concluding 
remarks are presented in Sec. IVII 



II. THE MODEL, THE WEIGHT DISTRIBUTIONS, ROTATED ELECTRONS, AND THE HOLON, 

SPINON, AND PSEUDOFERMION DESCRIPTION 

The ID Hubbard model reads, 

3, " 3 

where a and Cj jtT are spin-projection a =T,| electron operators at site j = 1,2, ...,N a and nj_ a — c] a c^ a - The 
model ([Ql describes N\ spin-up electrons and Ny spin-down electrons in a chain of N a sites. We use periodic boundary 
conditions and units such that the Planck constant and electronic lattice constant are one. We denote the electronic 
number by N = JVf + N± . The number of lattice sites N a is even and very large and in the units used here the chain 
length reads L = N a . We consider electronic densities n = n-\ + n| and spin densities m — n-\ — n\ in the range 
< n < 1 and < m < n, respectively, where n a — N a /L = N a /N a . In the thermodynamic limit the Fermi momenta 
read kp = unjl and kp a — Trn a . We denote the electronic charge by — e. 

We consider the following general fmite-w ./V-electron spectral-weight distribution, 

Bb(k,u>) = Y, \(f\6 l jv(k)\GS)\ 2 s(oj-l[E f -E GS }); lu> > ; I = ±1 . (2) 
/ 

Here the general A/"-electron operators Otf~(k) = 0\j-(k) and OTf 1 (k) = Oj^{k) carry momentum k, the / summation 
runs over the excited energy eigenstates, the energy Ef corresponds to these states, and Eqs is the ground-state 
energy. The local operator Oj^j = 0\j- j or Oj^j = Oj^.j is related to the corresponding momentum-representation 

operator O l j^{k) of Eq. by a Fourier transform. (Examples of such A/"-electron operators are given in Ref. (1C|-1 
For simplicity, we use in expression a momentum extended scheme such that k G (-co, +oo), yet it is a simple 
exercise to obtain the corresponding spectral function expressions for the first Brillouin zone. The general spectral- 
weight distribution (01 obeys the following sum rule, 
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dk dwB l M {k,u) = Y J {GS\dj} J 6 l u , j \GS); I = ±1 . (3) 
-oo JO „_j 



When j is hermitian one has that O 1 ^ ^ — Oj} j. In this case the convention is that I = +1 and thus the energy 
ui is positive. 

The concept of a rotated electron plays a key role in the pseudofermion description. Concerning its relation to the 
holons, spinons, and cO pseudofermions whose occupancy configurations describe the energy eigenstates of the model 
JU, see Ref. [Tl|. The holons have r;-spin 1/2, r^-spin projections ±1/2, charge ±2e, and spin zero. The spinons have 
spin 1/2, spin projections ±1/2, and no charge degrees of freedom. Throughout this paper we denote such holons 
and spinons according to their value of ry-spin projections ±1/2 and spin projections ±1/2, respectively. Furthermore, 
the cO pseudfermion is a composite quantum object of a chargeon of charge — e and an antichargeon of charge +e. It 
has no 7/-spin and no spin degrees of freedom. The 2^-holon (and 2z/-spinon) composite cv pseudofermions (and sv 
pseudofermions) are ?7-spin and spin zero (and spin zero) quantum objects (and with no charge degrees of freedom). 
These composite quantum objects contain an equal number v of —1/2 holons and +1/2 holons (and —1/2 spinons 
and +1/2 spinons). The ±1/2 holons (and ±1/2 spinons) which are not associated with 2z/-holon cv pseudofermions 
(and 2zy-spinon sv pseudofermions) are the ±1/2 Yang holons (and ±1/2 HL spinons). We denote the numbers of 
ctv pseudofermions and av pseudofermion holes by N av and N^ v , respectively, where a = c, s and v = 0, 1, 2, ... for 
a = c and v = 1,2, ... for a = s. (The value of N% v is given in Eqs. (B7) and (B8) of Ref. 0-) Moreover, L a ±i/ 2 
denotes the number of ±1/2 Yang holons (a = c) or ±1/2 HL spinons (a = s). 

Often in this paper we use the notation av ^ cO, si branches, which refers to all av branches except the cO and 
si branches. Moreover, the summations (and products) £ Ql/ , £ ai/=c0iSl! and £a^cO,«i (and Y[ av , Y[ av=cSS>al , and 
Ilai/^co si) run over all av branches with finite av pseudofermion occupancy in the corresponding state or subspace, 
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the cO and si branches only, and all av branches with finite av pseudofermion occupancy in the corresponding state 
or subspace except the cO and si branches, respectively. 

The pseudofermion description refers to a Hilbcrt subspace called the pseudofermion subspace (PS) in Rcf. |llj |. 
in which the A/"-electron excitations O l j^(k)\GS) are contained. The PS is spanned by the initial ground state and 
the excited energy eigenstates originated from it by creation, annihilation, and particle-hole processes involving the 
generation of a finite number of active pseudofermion scattering centers, —1/2 Yang holons, and —1/2 HL spinons plus 
a vanishing or small density of low- energy and small- momentum av = cO, si pseudofermion particle- hole processes. 
All these processes are defined in Refs. [10, EH 

The a^-pseudofermion discrete canonical-momentum values have a functional character and read, <jj — qj + 
QtAlj)/L = [2ir ; lL]If v + QtM])/L where j = l,2,...,N* v and N* av = N au + N% v . In this paper we use the 
designation functional for all quantities whose value depends on the energy eigenstates bare-momentum occupancy 
configurations. According to the studies of Ref. ^3], Qa vilj)/^ 1S a av pseudofermion scattering phase shift given 
by, 

N a<»' 

QLfe)/2 = ^E E *«*«^(<&><&0AJW(<&0; 3 = 1,2,~,K U . (4) 

ot'v' j' — l 

Here AN al/ (qj) = AAf av {qj) is the bare-momentum distribution function deviation AN au {qj) — N alJ {qj) — N® u (qj) 
corresponding to the excited energy eigenstate. (The ground-state densely-packed bare-momentum distribution func- 
tion N® u (qj) is defined in Eqs. (C.1)-(C.3) of Ref. This deviation is expressed in terms of the bare-momentum 
qj = [2ttI" v ]/L, which is carried by the av pseudoparticles, where I? u are the quantum numbers provided by the Bethe- 
ansatz solution Although the av pseudoparticles carry bare-momentum qj, one can also label the corresponding 
av pseudofermions by such a bare- momentum. When we refer to the pseudofermion bare- momentum qj , we mean that 
qj is the bare- momentum value that corresponds to the pseudofermion canonical momentum q~j = qj + Q® v (qj) / L. 
In contrast to the av pseudoparticles, the av pseudofermions have no residual-interaction energy terms. Instead, 
under the ground-state - excited-energy-eigenstate transitions the av pseudofermions and av pseudofermion holes 
undergo elementary scattering events with the a'v' pseudofermions and a'v' pseudofermion holes created in these 
transitions |12(. This leads to the elementary two-pseudofermion phase shifts n <fr al/ , Q v'l£ij l' j) on the right-hand side 
of the overall scattering phase shift I0J, defined by Eqs. (15) and (A1)-(A13) of Ref. Moreover, the overall av 

pseudofermion or hole phase shift, 

Q a ,(qj)/2 = Q a j2 + Qt(qj)/2, (5) 

plays an important role in the pseudofermion theory [12J. Here Qau{qj)/ L gives the shift in the discrete canonical- 
momentum value q~j that arises due to the transition from the ground state to an excited energy eigenstate and 
can have the values Q au /2 — 0, ±7r/2 and is defined in Eq. (11) of Ref. 1 lj . In this paper we use boundary 
conditions such that Q° av /2 = 0, — sgn(fc) 7r/2, where k is the excited-state momentum relative to that of the initial 
ground state. Here we assume that for the latter state A^/2 and N are odd and even numbers, respectively. Q° av /L 
give s the shift in the discrete bare-momentum value qj that arises as a result of the same transition. In reference 
|12| it is confirmed that Q al/ (qj)/2, Eq. (J5J, is a ais pseudofermion or hole overall phase shift associated with the 
ground-state - excited-energy-eigenstate transition. Its scattering part is the functional Q^ v (qj)/2 given in Eq. 10} . 
(In contrast to the pseudofermions, the —1/2 Yang holons and —1/2 HL spinons are scattering-less objects |l2j.) The 
av pseudofermion creation and annihilation operators fq. tCa/ and fq jjav , respectively, have exotic anticommutation 
relations, provided in Ref. These anticommutators involve the overall phase shifts © and play a key role in the 
spectral properties. 

For the ground state the pseudofermion, —1/2 Yang holon, and —1/2 HL spinon numbers are given by Y c o = N, 
N s i = N±, N alJ — i c , -1/2 — L s _x/2 — f° r av 7^ cOj si. We call A^ and the ground-state cO and si 
pseudofermion numbers, respectively. The ground-state av — cO, si bare-momentum distribution functions are 
such that there is pseudofermion occupancy for |g| < q% av and unoccupancy for q% av < \q\ < <7°„, where in the 
thermodynamic limit the Fermi-point values are given by, 

IfcO = 2k F ; QfsI = k Fi ■ (6) 
Moreover, for that state the limiting bare-momentum values of both the av — cO, si and av ^ cO, si bands read, 

Qco = * ; Isi = k F] ; <7c„ = [?r - 2fc F ] , v > ; q°„ = [k F1 - fc F J , v > 1 . (7) 
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Within the thermodynamic limit, when in this paper we refer to creation (or annihilation) of av = cO, si pseud- 
ofermions with bare-momentum values at the Fermi points, we mean that such a process leads to a densely packed 
bare-momentum pseudofermion distribution for the final state, where the created (or annihilated) pseudofermions oc- 
cupy neighboring discrete bare-momentum values ±q° FalJ + 2jn/L such that j — ±1, ±2, ... (or j — 0, =Fl, =f2, ...). Sim- 
ilarly, when we refer to creation of av ^ cO, si pseudofermions at the limiting bare-momentum values, we mean that 
such a process leads to a densely packed bare-momentum pseudofermion hole distribution for the final state, where the 
created pseudofermions occupy neighboring discrete bare-momentum values ±q^ v + 2jir/L such that j = 0, =pl, =p2, .... 
As further discussed below, the cv pseudofermions of limiting bare-momentum value q = ±q° v = ±{ir — 2kp] and 
sv pseudofermions of limiting bare-momentum value q = ±c/j? = ±[^ft — &fjJ can be described in terms of 2v 
independent holons and 2v independent spinons, respectively |l2j . 

Our study of the metallic phase of the model Q involves several numbers and number deviations. N^ NF is 
the number of finite-momentum and finite-energy av = cO, si pseudofermion particle-hole processes. The quantum 
number i = sgn(g)l = ±1 refers to the right pseudfermion movers [i = +1) and left pseudfermion movers (t = — 1) 
and AN F V L such that AN F U ±i is the deviation in the number of av pseudofermions at the right (+1) and left (— 1) 
Fermi points. In turn, the deviation in the number of av = cO, si pseudofermions created or annihilated away from 
these points is denoted by AN^f . The actual number of av pseudofermions created or annihilated at the right (+1) 
and left (-1) Fermi points is denoted by AN°^ ±1 . It is such that AJV^ L = AN^ F L + 1 Q° a „/2ir, where Q° av /2 is the 
scattering- less phase shift on the right-hand side of Eq. (JSJ. Furthermore, N F V t refers to the av ^ cO, si branches 
and is the number of av pseudofermions of limiting bare momentum q = Lq® ail such that t = ±1. The number of av 
pseudofermions created away from the limiting bare- momentum values is called N^ F . 



III. SPECTRAL- WEIGHT ELEMENTARY PROCESSES FOR DENSITIES < n< 1 AND < m < n 

In this section we study the pseudofermion occupancy configurations of the PS energy eigenstates. As in reference 
|llj |. we classify the pseudofermion elementary processes that generate the PS from the ground state into three types: 

(A) Finite-energy and finite- momentum elementary cO and si pseudofermion processes away from the corresponding 
Fermi points involving creation or annihilation of a finite number of pseudofermions plus creation of av ^ cO, si 
pseudofermions with bare-momentum values different from the limiting bare-momentum values ±9^; 

(B) Zero-energy and finite-momentum processes that change the number of cO and si pseudofermions at the 
i = sgn(g) 1 = +1 right and i = sgn(g) 1 = —1 left cO and si Fermi points - these processes transform the ground- 
state densely packed bare-momentum occupancy configuration into an excited-state densely packed bare-momentum 
occupancy configuration. Furthermore, creation of a finite number of independent —1/2 holons and independent 
— 1/2 spinons, including —1/2 Yang holons, —1/2 HL spinons, and —1/2 holons and —1/2 spinons associated with cv 
pseudofermions of limiting bare momentum q = ±g^ = ±[7r— 2kp ] and sv pseudofermions of limiting bare momentum 
q = ±<?g„ = ±[fc_F| ~ k-Fj.], respectively; 

(C) Low-energy and small-momentum elementary cO and si pseudofermion particle-hole processes in the vicinity of 
the l = sgn(q) 1 = +1 right and i = sgn(q) 1 = —1 left cO and si Fermi points, relative to the excited-state av = cO, si 
pseudofermion densely packed bare-momentum occupancy configurations generated by the above elementary processes 
(B). 

It is found in reference that the invariance under the electron - rotated-electron unitary transformation of the 
av pseudofermions created at limiting bare momentum q = ±<7°„ and belonging to av ^ cO, si branches implies that 
each of such cv pseudofermions (and sv pseudofermions) separates into 2v independent holons (and 2v independent 
spinons) and a cv (and sv) FP scattering center. By independent holons and spinons it is meant those which remain 
invariant under the electron - rotated-electron unitary transformation. (The Yang holons and HL spinons are also 
independent holons and spinons, respectively.) The above designation FP stands for Fermi points. Indeed, it is found 
in the same reference that the cO and si pseudofermion and pseudofermion hole scatterers feel the created cv (and sv) 
FP scattering centers as being cO (and cO and si) pseudofermion scattering centers at the Fermi points. Therefore, 
creation oiav pseudofermions of limiting bare momentum q — ±q& v and belonging to av ^ cO, si branches is included 
above in the elementary processes (B). 

In addition to the PS, the following subspaces play an important role in our studies: 

A CPHS ensemble subspace is a subspace of a canonical ensemble subspace with fixed and N± electronic numbers, 
which is spanned by all energy eigenstates with fixed values for the —1/2 Yang holon number L C) _i/2, —1/2 HL spinon 



■5 



number L s _ _i/2, cO pscudofermion number N c0 , and for the sets of av ^ cO pscudofcrmion numbers {N clJ } and {N sv } 
corresponding to the u =1,2, ... branches. Here CPHS stands for cO pscudofermion, holon, and spinon; 

A J- CPHS ensemble subspace is a subspace of a CPHS ensemble subspace spanned by the excited energy eigenstates 
with the same values for the numbers N^ Q NF , Nf 1 NF , AN F 0+1 , AN F lt AN^ +1 , ANf x _ l7 and sets of numbers 
{N F V +1 } and {N F V _ x } for the av 7^ cO, si branches with finite pscudofermion occupancy in the CPHS ensemble 
subspace; 

A J-CPHS subspace is a J-CPHS ensemble subspace whose numbers N^ NF and N^ NF have finite values; 

A reduced J-CPHS subspace is a subspace of a J-CPHS subspace spanned by all the energy eigenstates of the latter 
subspace which are generated from the ground state by elementary processes (A) and (B) only. 

The CPHS ensemble subspace dimension is given by, 



— - f3 ( N ° - -'?•>■ w -i n S n O • 



' x ' u>0 x ' i/'>1 v 

It is a product of the number of occupancy configurations of the cO pseudofcrmions, si pseudofcrmions, and av 7^ cO, si 
pseudofermions belonging to branches with CPHS ensemble subspace finite occupancy. Once the ground-state numbers 
Af? = N° and = N® are known and the corresponding excited-state values of these numbers are given by 
N c o = N® Q + AN c q and N s \ = + AN S \, one can define a CPHS ensemble subspace by the set of deviation number 
values AN c0 and AN s i and set of number values L Ct -1/2, L s , -1/2, an d {N al/ } for the av 7^ cO, si branches. 

Since the PS excited energy eigenstates involve the creation of none or of a finite number of av pseudofermions 
belonging to the av 7^ cO, si branches such that AN av — N av , all the CPHS ensemble subspace occupancy 

configurations of each of the latter av bands correspond to PS excited energy eigenstates. The same applies to the 
CPHS ensemble subspace cO (and si) pseudofermion occupancy configurations provided that n = 1 (and m = 0) 
for the initial ground state. Indeed, for such a density there are no cO (and si) pseudofermion holes for that state, 
and again all excited-state CPHS ensemble subspace configurations can be reached by a finite number of cO (and si) 
pseudofermion processes. However, for densities < n < 1 and < m < n the ground state has both "particles" and 
"holes" in the cO and si pseudofermion bands, whose numbers are such that N av — > 00 and N% v — > 00 as N a — ► 00. 
Thus, reaching from the ground state many of the excited-state CPHS ensemble subspace cO and si pseudofermion 
occupancy configurations involves in the thermodynamic limit an infinite number of finite-momentum and finite-energy 
pseudofermion processes. It follows that for the metallic phase at finite spin density only the J-CPHS subspaces of 
each CPHS ensemble subspace are contained in the PS. Indeed, J-CPHS ensemble subspaces with N^ NF and/or 
j^phNF mnmtc ^ j not belong to the PS and have no finite overlap with the one- and two-electron excitations. 



A. GENERAL DISTRIBUTION-FUNCTION DEVIATIONS FOR DENSITIES < n< 1 AND < m < n 

A J-CPHS subspace is spanned by the excited energy eigenstates contained in the corresponding reduced J-CPHS 
subspace and by all excited energy eigenstates generated from the former states by the elementary processes (C). Here 
we study the pseudofermion occupancy configurations of the excited energy eigenstates that span a J-CPHS subspace. 
There is a particular type of reduced J-CPHS subspace which we call point- sub space. A point-subspace is spanned by 
a single excited energy eigenstate whose deviation numbers and numbers are such that, 



AN au = AN F V = ]T AN F Vt L , av = cO, si ; N av = N F V = ^ N F ^ t , av ± cO, si . (9) 
t =±i i.=±i 

The number deviation AN F V obeys the relation AN au — AN^ F + AN F V , where AJV^/ and AN F V refer to the cO and 
si pseudofermion number deviations generated by the elementary processes (A) and (B), respectively. For the av = 
cO, si branches, the deviation number AN F V L is related to the current number deviation, AJ F V = | X^iiW ^-N F V L 
such that AN F v l — AN F v /2 + lAJ f v . Furthermore, for the av 7^ cO, si branches the number N F V L is associated 

with the current number J F V — 5 J2i=±i( L ) ^av, 1 sucn tnat ^av, 1 = N F v /2 + 1 J F V . 

The excited energy eigenstates that span a reduced J-CPHS subspace have well-defined cO and si pseudofermion 
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right (i = +1) and left (t = —1) bare-momentum and canonical-momentum Fermi points, 



lFau,i = ^F av + ^QFau,i; ^QFau, l = 4 y- AiV^ , ; OJ/ = cO, si ; t = ±1 , (10) 

and 



gFa„, t = ^F Q , + Ag Fai/ , t ; Ag Fcw , t = t y AiV^ , + Q ^ L ^^ ■ a v = cO,sl; t = ±l, (11) 

respectively, where q FQ ,„ is the initial ground-state Fermi-point value given in Eq. JBJ. Here Q*„(g)/2 is the cw 
pseudofermion overall scattering phase shift (@J. We call 2A L aiy the square of the av, i Fermi-point value shift in units 
of2n/L, 

2A^ ee (^f) 2 - + 9sMs^ = ^AJC, t + g^M^l ) 2 ; W = c0,sl; t = ±l. (12) 

This quantity plays a key-role in the spectral properties, as discussed in Sees. IV and V. It has the form given in Eq. 
Ill 21 1 for A/"-electron spectral functions leading to deviations in the electronic numbers, A7V, such that |AJV| = 0, 1, 2. 
Our analysis is limited to such spectral functions. Our studies of these functions involve the use of specific expressions 
for the pseudofermion bare-momentum distribution function deviations. The cO and si pseudofermion bare-momentum 
distribution-function deviations of the excited energy eigenstates that span a J-CPHS subspace have the following 
general form, 



AN ap {q) = AN^ F (q) + AN F v {q) + AN^; F (q) ; av = cO, si . (13) 

The bare-momentum distribution- function deviation AN^f (q) corresponds to the above elementary processes (A) 
and is given by, 



AN^(q)= S gn(AN^) 



2 k 



E 

i=l 



6(q 



2tt 

T 



[S(q- qj )-d(q-q' 



i=i 



av 



cO, si . 



(14) 



where sgn(AN^f) = sgn(A7V at/ ) and q t with i = 1, |AiV^/| gives the set of bare-momentum values associated with 
creation (AN^ F > 0) or annihilation (AN^ F < 0) of av pseudofermions away from the Fermi points. Moreover, qj 
and q'j with j = 1, N^l NF denote the bare- momentum values associated with the particles and holes, respectively, 
of the finite-energy and finite-momentum av pseudofermion particle-hole processes. 

The bare-momentum distribution function deviation AN F v (q) results from the above elementary processes (B) and 
reads, 



A<^) = T <KM-4 Q JAAC,Ai,sgn (9) i; av = cO,sl. (15) 

Finally, the bare- momentum distribution function deviation AN%,„ (q) is associated with the elementary processes 
(C) and is given by, 



AN P 



h F 



(?) 



2tt 



N phF 

E 

(.=±i j L =i 



av = cO, si , 



sgn(g) 1 . 



(16) 



In the latter equation qj L and * with j L — 1, N^ F denote the bare- momentum values associated with the 
" particles" and " holes" , respectively, in the vicinity of the Fermi point of Eq. IjlOl) corresponding to the av, i sub- 
branch index value i = sgn(g) 1 = ±1. Such bare-momentum values are ordered as uq\ l < LQ2 l < ... < iq NP hF and 



< i -i , 2 l < ■■■ < L i' 



a'; 



for av — cO, si and L = ±1. Once the elementary processes (C) generate the J-CPHS 
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subspaces from the corresponding reduced subspaces, the set of "particle" and "hole" bare-momentum values {qj L } 
and {q'j L } are such that q = qj t and q = q' ^ was unoccupied and occupied in the corresponding initial excited energy 
eigenstate belonging to the latter subspace, respectively. The number N^ F is such that N%£ F < m av ,L, where m av , t 
is the number of av = cO, si pseudofermion elementary particle-hole processes of momentum l[2tt/L] involved in the 
deviation AN££ F (q) of Eq. (jlfijl ■ Different such deviations may correspond to the same number wn av , l ■ 

The av ^ cO, si branches have no finite pseudofermion occupancy for the ground state. The av ^ cO, si pseud- 
ofermion occupancy configurations of the excited energy eigenstates which span a J-CPHS subspace can be defined by 
the values of the sets of numbers {N F V L } associated with the independent holons or spinons and the av FP scattering 
centers and the following bare-momentum distribution-function deviations, 

^f(q) = i q' z ) ; av ± cO, si . (17) 

t=i 

Here q\, with i = 1, JV^f, gives the set of bare-momentum values of the pseudofermions created by the elementary 
processes (A) away from the limiting bare- momentum values ±g°„. 

B. THE ENERGY AND MOMENTUM SPECTRA FOR DENSITIES < n < 1 AND < m < n 

Our studies of Sees. IV and V require an analysis of the energy and momentum spectra beyond that of Ref. [Io| . 
The pseudofermion dispersions e au (q) determine the form of the energy excitation s pec trum of the weight distributions 
studied in these sections. These energy bands are plotted in Figs. 6-9 of Ref. [l7| as a function of q for m — > 
and several values of U/t and n. By suitable manipulation of Eqs. (C15)-(C21) of Ref. 0, we find the following 
dependence on the pseudofermion bare-momentum q, 

U r +Q ~ 

e c0 (g) = M ~ Mo# — — 2t cos k°(q) + 2t dk $ c0 , c0 (k, k°(q)) sin k , (18) 
1 J-Q 



, (q) = 2/iv + e° cv (q) = {2 l i-U)v + At Re { ^1 - (a°,(<z) - iz/^) 2 } + 2t | ^ dfc <J> c0 , cv (k, A° cv (q)) sin k, (19) 



and 



e sv (q)=2vfioH + e sv (q)=2vfx H + 2t dk $ c o,sv(k, A° v (g)) sin k . (20) 

Here the two-pseudofermion phase shifts $ are given in Eqs. (|A3|) and l|A4|l of Appendix A, Q is the parameter 
defined by Eq. 1A5(I of that Appendix, and fj, — fi(n, to, U/t) and H = H(n, m, U /t) can be expressed in terms of the 
same phase shifts as, 



u r +Q ~ r +Q ~ 

fj. = tt + 2t cos Q - t / dk [2$ c0 , c0 (k,Q) + $ c0 , „i (k, B)\ sin & ; fi H = -t / dfc $ c0 . s i (fc, B) sin fc . (21) 

* J-Q J-Q 



The ground-state rapidity functions k°(q), A® u (q), and A°„(g) appearing on the right-hand side of Eqs. l(T5|l - lf2l7|) 
are defined in terms of their inverse functions in Eqs. (|A1I) - (|A2|) of Appendix A. (Analytical expressions for the 
rapidity functions k°(q), A°„(g), and A° sl/ (q) as to — > and both U/t — > and f//i >> 1 are provided in Ref. [lq.') 
The zero-energy, 2/w-energy, and 2v/j 1 c,H-energy levels of the pseudofermion energy bands (|18fl - H2Ufl relative to the 
ground-state energy are such that, 

e cQ (±2k F ) — e s i(±kpi) = 0; e cv (±[iT - 2k F ]) = 2vfi , v>0; e sv {±[kF] - fcj?;]) = 2vfx H , v>l: 

e° w (±[7r - 2k F }) = 0, ^>0; £,(±[4^ - = , i/ > 1 . (22) 
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Creation onto the ground state of one cv pseudofermion or sv pseudofermion belonging to v > and v > 1 
branches, respectively, is a finite-energy process whose minimal energy value is given by e cl/ (±[7r — 2kp]) = 2v/i or 
£si/(0) = 2i>[i H + e^(0), respectively. Creation onto the ground state of one independent —1/2 holon or independent 
— 1/2 spinon requires energy 2/i or 2/io-ff, respectively 9]. In the limit n — *■ 1 (and m — > 0) both the energy and 
bare-momentum band widths of the dispersions e C v(q) for v > (and e s „(q) for v > 1) vanish. Moreover, while the 
energy band width At and bare- momentum band width 27r of the dispersion e c o(q) remain the same for all values of 
U/t, the energy band width of the dispersions of all other pseudofermion branches vanishes as U/t — > oo. This effect 
follows from the localized character reached by the av pseudofermions of v = 1,2, ... branches as U/t — > oo 0,^3- 
In that limit the ^ cO pseudofermions and sz^ pseudofermions separate into 2v independent holons and spinons, 
respectively. 

The group velocity v av (q) and the light velocity v av are given by, 

Vau{i) = ^"q^ > a11 branches; v au = v av (q° Fav ) , av = cO, si , (23) 

and play an important role in the spectral-weight distribution expressions obtained in Sec. V. 

The cw-branch finite-energy spectrum associated with the elementary processes (A) corresponds to the reduced 
J-CPHS subspaces and reads, 



\ olu \ a is 



av = cO, si : 

i=l j=i 

AE^dq't}) = e^iq'i) ; m/ ^ cO, si , (24) 



where the pseudofermion energy bands are given in Eqs. (JTSJ-lj^DJ. 

The finite-momentum spectrum has contributions both from the elementary processes (A) and (B) and also corre- 
sponds to the reduced J-CPHS subspaces. It is given by, 

AP av ({ qi }, { qj }, tf.}) = AP£ + AP?f (M, { Qj }, tfj}) ; av = cO, si ; 

*Pc*(Wi}) = k w + AP c N /({ q ' t }); cv^cQ; AP st ,({g'J) = £ q\ ; sv + s\ . (25) 

i=i 

Here AP£ V is given in Eq. (67) of Ref. 10], 



AP£f({ qi },{ qj }, {q'A) 



sgn(AAr) 



\AN NF \ 

E 

i=l 



N phNF 



3j - 5 



cO, si . 



(26) 



AP c ^({g'J) = fc c ,= (l + ^Ar c ^; cv ^ cO , (27) 
i=l 

and g^Q,^ is the cO and si pseudofermion ground-state Fermi-point value given in Eq. ©. The general finite-energy 
and finite-momentum spectrum generated by the elementary processes (A) and (B) reads, 



AE = AE({q i },{q' i },{q j },{q' j }) = £ AE av ({qi}, {qj}, {q'j}) 



af—cO. si 



E 

av^cO, si 



E 

Ct—C. s 



(28) 



and 
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AP = AP({q i },{q' i },{q j },{ q ' j }) = P + £ AP^/ ({%},{<?,}, {<?',}) 

CM/— cO, Si 

+ £ AP c f({l'J)+ E Ap -({^})» 

ci/^cO sv^sl 



(29) 



where 



P = 



Pa 



cv^cO 



if* - 
«, — 



AP F 



(30) 



af=cO, si 



correspond to the independent holons (a = 
J, respectively. Both the spectra (|28|l and 



= c) and spinons (a = c) and are 
l(29"|) are associated with reduced 



and the energy E a and momentum P a 
given in Eqs. (66) and (67) of Ref. ^ 
J-CPHS subspaces. 

Generation of the J-CPHS subspaces from the corresponding reduced subspaces involves the elementary processes 
(C). Given the linear av = cO, si pseudofermion energy dispersion near the Fermi points, these processes lead to 
small momentum and energy values such that, 



k' 



E L T m ' c 

au—cO, si t=±l 



2tt 



aiy—cO, si l—±1 



(31) 



Here m av , i is the number of elementary av pseudofermion particle-hole processes of momentum l[2tt/L] defined above 
and v av is the av = cO, si pseudofermion light velocity given in Eq. (|23(l . Thus, the elementary processes (C) generate 
a set of excited energy eigenstates with energy and momentum given by those of the initial reduced-subspace state, 
Eqs. 1)281) and (|29|) . respectively, plus the small energy and momentum given in Eq. H31|) . According to the results 
of Ref. the elementary processes (C) have a non-interacting character in terms of av = cO, si pseudoparticles. 
In turn, for the remaining low-energy excitations such pseudoparticles have residual interactions, in contrast to the 
corresponding av — cO, si pseudofermions. This pseudoparticle non- interacting character implies that the energy 
spectrum of Eq. (|31|l remains linear in m av _ L for small finite values of to q „. L /N a as N a — > oo. 

The energy ujq for transitions from the ground state to a point-subspace plays an important role in the spectral- 
weight distributions. Such an energy is that of the excited energy eigenstate whose deviation numbers are given in 
Eq. © and reads, 



LOq 



E E a = 2 M M C; _ 1/2 + 2fi H [M s> _ 1/2 - N sl ] 



(32) 



where N s % = + AN s i. We note that each CPHS ensemble subspace contains at least one point-subspace and that 
all point-subspaces of a CPHS ensemble subspace have the same energy, given in Eq. I|32l) . Furthermore, the smallest 
excitation energy for transitions from the ground state to a CPHS ensemble subspace reads, 



^CPHS 



sv^sl 



(33) 



where ujq is the energy l|32() of the point-subspace(s) contained in the CPHS ensemble subspace and e%(0) < for 
77i > and e°„(0) — » as m — > 0. Thus, for m — > 0, one has that E^ PHS = ujo = 2/iM c _ 1 / 2 - In this case the 
energy value E^ PHS = ojq — corresponds to a M c _ 1 ^ 2 = CPHS ensemble subspace, E^ PHS = ojq = 2/x to the 
M c> _i/2 = 1 first-Hubbard band CPHS ensemble subspace, Eq PHS = ujq = 4/x to the M C) _i/2 = 2 second- Hubbard 
band CPHS ensemble subspace, and so on. We recall that the number of —1/2 holons, M c , -i/2 ; equals that of 
rotated-electron doubly occupied sites 

The excited energy eigenstate that spans a point-subspace has energy u = Iujq and momentum k = Ikn relative to 
the initial ground state, where kg is given by, 



ko 



E 



Pa = n[L 



c, -1/2 



4k f 



2fe 



ajs 



2E J - 



(34) 
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Here k F is provided in Eq. (|3U[) . A property of the one- and two-electron spectral- weight distributions is that for 
small values of the energy deviation (u) — Iujq) and m — > there is spectral weight only for values of momentum k such 
that (k — lko) is also small. While the momentum spectrum tt [L Cj -1/2 + S^Li v -^cu\ i s additive in the momentum 7r of 



refer 



each independent —1/2 holon, the momenta Akp ^J^ + YluLi Jcv + J2T=2 ^su an d 2&fj.[ AJ,^ — 2^^L 2 J* 

to the cO and si pseudofermion Fermi points current number deviations. Note that although the current numbers J^ v 
and J\ v arise from the cv 7^ cO and av 7^ si FP scattering centers, they lead to momentum contributions associated 
with the cO Fermi points and cO and si Fermi points, respectively. These current numbers correspond to phase 
shifts whose scattering centers are felt by the cO and si pseudofermions has being at the av = cO, si Fermi points 
|l2j| . Point-subspaces contained in the same CPHS ensemble subspace have the same energy but different momentum 
values. For a given CPHS ensemble subspace, we call first point-subspaces and second point-subspaces the subspaces 
with the momentum l)34|) given by k = ±|fc | which correspond to the minimum value of |fcol an d the second smallest 
value, respectively. 

C. THE ELEMENTARY PROCESSES AND THE ACTIVE SCATTERING CENTERS 

Active scattering centers are those which contribute to the scattering phase shift (0J. We find below that only the 
elementary processes (A) and (B) generate active scattering centers. For the excited energy eigenstates described 
by the bare- momentum distribution function deviations given in Eqs. (|13fl -i|17 [) the scattering phase shift Q can be 
written as, 

QL(q)/2 = Qtl NF) (q)/2 + Qtl F) (q)/2 ■ (35) 

Here the phase shifts Q%i NF \q)/2 and Q%i {q)/2 result from active scattering centers generated by the elementary 
processes (A) and (B), respectively. The spectral properties are mainly controlled by the scattering phase shifts 
QavitQFau)/^ °f the av = cO, si pseudofermion and hole scatterers of bare momentum values q « ±q Fav . The 
scattering phase shift Q%i NF \iq Fav ) /2, where av = cO, si and i = ±1, can be written as, 

Qti NF) (ui Fau ,{<ii}d<i'i},{<ij},Wj}) r t fA MNF^sr?' 1 * , o n 

a' y' — cO, si i— 1 

pfPhNF N N , F , 

+ [®<*>>,0SsfaFav>Q3)-$«v, a >u>(t4F a vi</j)j+ ® a »' a ' v ' ^Fau > l' i) } • ( 36 ) 

In turn, the phase shift Qai {tq^ av )/2 was derived in Ref. and reads, 



AN F AN F 

^ iv c0 | H) ^ v sl 
1 Sqi/ si 



AA' 



2 -'^^ 2 

oo 

Cav si AJsl ~ ^ J sv 



AJav + Cay cO ^^cO + ^cv + ^ 



(37) 



Here the Fermi-point parameters CawaV are symmetrical (j — 0) and antisymmetrical (j = 1) combinations in 
the a'v' — cO, si pseudofermion scattering centers right and left Fermi-point values of the elementary av = cO, si 
two-pseudofermion phase shifts in units of 7r given by, 



t.aua'v> = S <*,<*' + X! ^^^^^'^'^Fav^lFa'v') 5 .7=0,1; CiV = cO, si ; aV = cO, si . (38) 

t=±l 

The parameters l|38|l appear in the spectral- function exponents derived in Sees. IV and V. In the limit m —* 0, 
they are given by £? 0c0 = l/f , Crfisi = 0, Csico = -1/V2, Csisi = V^, CcOcO = Co, CcOsi = Co/2, ^ lc0 = 0, and 
Csisi — 1/^/2- Here £o is the parameter defined in Eq. (74) of Ref. 01 an d in the text above that equation. It is 
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such that £o ~^ V% and £o - * 1 as U/t — ► and U/t — * oo, respectively. As further discussed in Ref. [l^, the form of 
the general phase-shift expression l|37|) confirms that although the current numbers J F V and Jf v are associated with 
the cv ^ cO and sv ^ si FP scattering centers, the cO and si pseudofermions feel such centers as being at the cO 
Fermi points and cO and si Fermi points, respectively. 

The following functional plays a major role in the spectral properties, 

Co = Co(te},{</J, {?,•}, {?',•}) =2A c0 + 2A sl ; 2A a „ = 2A+3 + 2A"* ; au = cO,sl. (39) 

Moreover, by use of of the general expressions i|3T)|) and (|3TJl one finds that the four functionals 2A L ain Eq. lfE2)l. can 
be written as follows, 



2A*„ = 2A 



T ''Scti/eO 



2tt 



A /V F A N F 



co [ A Jfo + E J £ + E J -l + & »i [ A J i"i - 2 E J -l ) i = cO, si , t = ±1 , (40) 



i/=2 



!/=2 



where the phase shift Q a l ( ic lFav) is that of Eq. 

Analysis of the above general phase-shift expressions confirms that the cO and si elementary processes (C) associated 
with the bare-momentum distribution function deviations ANP^ F (q) of Eq. (|16(1 do not generate active scattering 
centers. Let us consider a "hole" and "particle" created at bare-momentum values q' and q' + [2n/L]Nf, respectively, 
where Nf — tl, l2, ... is a finite integer number and i = ±1. According to Eq. QJ, such a elementary process leads 
to the scattering canonical-momentum shift Q% v {q)/L = [2n/L] [$ av , a'v'il, <i + [2^/L]Nf) — ^ a v,a'u'(<l,Q J )] such 
that Q%„(q)/L w [2n/L] 2 Nf [d^ a v,a'v'(QiQ')/dq'], for all pseudofermions and holes whose branches have finite 
pseudofermion occupancy in the excited energy eigenstate. Here the second expression corresponds to the leading- 
order contribution in 1/ L. Since for the pseudofermion description discrete bare-momentum and canonical-momentum 
contributions of order [l/£p such that j > 1 have no physical significance we thus conclude that Q^ v (q)/L = 0. 

Once all the N^ F "particle" and "hole" pairs of the bare- momentum distribution function deviation AN^ F (q) of 
Eq. H16|) are such that qj L — q' a = [2tt /L]Nf where Nf — tl, t2, ... is a finite integer number, such a deviation does not 
contribute to the scattering phase shifts. Therefore, up to first order in 1/L, the scattering phase-shift contributions 
from the deviations involving the summations of the N%£ L "particle" S- functions [2ir/L] 5(q — qj t ) of the deviation 
(|16[1 are exactly canceled by those involving the N^ F t "hole" 5-functions — [2tv/L] S(q — q'j t ). That the av = cO, si 
pseudofermion "particle" and "hole" scattering centers generated by the elementary processes (C) are not active 
scattering centers implies the validity of the following property: The overall scattering phase shift overall phase 
shift ©, and functional l|4U|) have the same value for the whole set of J-CPHS subspace energy eigenstates generated 
by the elementary processes (C) from a given energy eigenstate of the corresponding reduced J-CPHS subspace. This 
property plays an important role in our study. 



IV. FINITE-ENERGY WEIGHT DISTRIBUTIONS FOR DENSITIES < n < 1 AND < m < n 

Our starting point for the derivation of closed-form analytical expressions for the one- and two-electron weight 
distributions obtained in the ensuing section is the general A/"-electron spectral function given in Eq. (36) of Ref. [Ig, 



i=0 



E 



B l '\k, lu) 



= 1 



l = ±l. 



(41) 



where for the present case of densities < n < 1 and < m < n the summation over the numbers NP^ NF is limited to 
finite values and thus refers to reduced J-CPHS subspaces only. On the right-hand side of Eq. 141(1. c\ is the constant 
of the operator expressions given in Eqs. (32)-(34) of Ref. [Hj such that c\ — > as U/t — > co for i > and the 
function B l,l (k, u>) is defined in Eq. (68) of that reference. Let us rewrite the av pseudofermion spectral functions 
introduced in Eq. (63) of Ref. [lj| as follows, 
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B^(k,u) = N a {(e(AN»f) n [^( E + E )]+e(-AJC') II E 

x II [^r( E + E )]v" E ]}*(w-^^)**,iap ,; ai/ = cO,*l; 
B^i^u) = N a [ f[ — ]s{u-l&E av )6 k<lAPa „; av^cO,sl; l = ±l; t = 0,1,2,..., (42) 



1=1 



where the energy and momentum spectra are given in Eqs. (|24|l and H25|) . respectively. On the right-hand side of 
Eq. I|42l) and in all expressions given below the 8 function is such that Q(x) = for x < and 0(x) = 1 for x > 0. 
Moreover, here we have explicitly written the summation Ylj-cPHS-av-(A) °f expression (63) of Ref. E3, which 
runs over the reduced J-CPHS subspace av pseudofermion occupancy configurations generated by the elementary 
processes (A). The number of such occupancy configurations is given by Y\ av D au , where for the av = cO, si branches 
the dimension D av reads, 



D av = {e(sgn(A</)) ( A ^ F ) + e(-sgn(A</)) } (^L) (^,) , (43) 

and is given in Eq. (64) of Ref. [13 for the av ^ cO, si branches. We recall that Eq. (|43|l refers to densities < n < 1 
and < m < n and thus the numbers N c0 , N^ , N sl , and are such that iV c o — > oo, N^ Q — > oo, N s i — > oo, and 
Njji — > oo as 7V a — > oo, whereas for the PS excited energy eigenstates the numbers AN^ U F , AN F V , and N^ NF have 
finite values. Thus, in expression (|43|) we have not considered small finite corrections to the numbers N c o, N^ , N s i, 
and Ng 1; which lead to higher order corrections that vanish in the thermodynamic limit. Moreover, in expression 
(14211 wc could treat the bare-momentum summations as independent for each created or annihilated pseudofermion. 
It follows that for N a — > oo, exclusion of the occupancy configurations where two or several created or annihilated 
pseudofermions meet at the same discrete bare-momentum value leads to vanishing corrections to expression (|42H . 

From use of the spectral function expression (|42|) in the second expression of Eq. (68) of Ref. ^(| > we arrive to the 
following expression for the function B l ^(k, us) appearing on the right-hand side of Eq. 14111 . 



n [(e(AA&*) n lw( e + e )]+e ( -A</) n lw e 

- ' " 1 ' ' i=1 a q^=-q F , 

e(n-i[w-iAE]\ 



ai^=c0 1 si 



N phNF 



n ht( e + e )£ e ]]}[ n e 

l' j = -1 Fav av^c0,sl i=l q' i = —q[ 



e(z[w-zAP]) e 



Z[w - ZAP] 
life- ZAPI 



u a i ) B M ( w - ZAP 



w-ZAP\ 
k - ZAP J 



i = 0,1,2,..., Z = ±l. 



(44) 



Here the general energy functional AP = AE({qi}, {q'i}, Wj}) ancl momentum functional AP = 
AP({<7j}, {q'i}, {qj}, W j)) are given in Eqs. (|2*S|) and lt2$|) . respectively, and in the thermodynamic limit the value 
of the small positive energy Q, which corresponds to the energy range of the elementary processes (C), is controlled 
by the functional £o, Eq. (|39() . of the corresponding initial rcduced-subspace state. (The energy £1 vanishes when 
Co — > and the weight distribution becomes (5-function like, as further discussed in Sec. V.) In numerical studies of 
the spectral-function expressions derived in this paper the small energy f2 is treated as a parameter whose value for 
fixed finite values of n, m, and U/t is determined by imposing the exact sum rule J3J in the end of the calculations. 
The bare-momentum summations of expression (|44() run over all Y\ au ^ aL ' 0CCU Pancy configurations generated by the 
elementary processes (A). For the whole (k, w)-plane and except for k and u values such that uj ~ i v av (k — Iky) + Iujq, 
where av — cO, si, b = ±1, ujo is provided in Eq. I|32l) . and ko in Eq. (|34ll . the function B l ' l {Au>,v) on right-hand 
side of Eq. 1441) is the convolution of the cO and si spectral functions given in Eq. (62) of Ref. These functions 
read, 
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pM(Au,, v) = ± J du>J2 B% sl (M« - k', Au - J) B l l a (V, <J) 

00 k' 

'' 1 ' dJ / dk ' B Q sl ( A ^/ V ~ Aw - W 'J 5 Q c0 ( fc '' W ) ' 



•/ — sgn(v) Aa;/i? c o 

B$ au {k',u,') = ^ KO|^-G*a.^p-h, w F! G s,^ (45) 

J-CPHS-uv~(C) 

In this equation i = 0, 1, 2, ...,/ = ±1, Acj = (w — IAE), Ak = (k— ZAP), the summation ^ j-c'PHS-av-(c) runs over 
the J-CPHS ensemble subspace av = cO, si pseudofermion occupancy configurations generated by the elementary 
processes (C), the matrix-element generators are given in Appendix B, and. vn^u = J2l=±i m au, i- The identity of the 
first and second expressions of Eq. I|45[l follows from the structure of the spectrum of Eq. <|31(l . which is generated by 
the elementary processes (C). However, we note that the validity of the second expression and of the spectral- function 
expressions given below is limited to electronic densities such that v c q > v s \. We emphasize that this inequality holds 
for all electronic densities n of the metallic phase except for a small domain in the vicinity of n = 1. The velocity v 
appearing in the argument of the function l|45|l plays an important role in our study and is given by, 



Aw (u - IAE) . . „ / a j \ 7 

v= A~k = Fib - zap) ' Sgn ^ = s s n ( A/c ) 1 ; M > v si ■ ( 46 ) 

The inequality \v\ > v s \ and the theta- functions appearing on the right-hand side of Eq. 144|) follow again from the 
structure of the spectrum of Eq. (f!?T|) . 

Each excited energy eigenstate generated from the ground state by the elementary processes (A) and (B) corresponds 
to one point, (ZAP, IAE), of the finite-weight (k, w)-plane region associated with such a subspace. The set of all such 
points generated by the bare-momentum summations on the right-hand side of Eq. I|44(l span a well-defined finite- 
weight (k, cj)-plane point-like, line-like, or surface-like domain. From the reduced J-CPHS subspace excited energy 
eigenstates corresponding to points (ZAP, ZAP) in the vicinity of the (k, w)-plane point that the spectral-function 
expression (I44|l refers to, the elementary processes (C) generate excited energy eigenstates whose momentum and 
energy relative to the initial ground state is k and u), respectively. Thus, the momentum ZAP and energy ZAP of the 
reduced J-CPHS subspace excited energy eigenstates relative to the initial ground state are such that (u> — IAE) and 
(k — ZAP) are small and \v\ = \(u> — lAE)/(k — ZAP)| obeys the inequality \v\ > v s \. 

The av = cO, si pseudofermion spectral function B 1 q (k' , u>') of Eq. I|45() can be expressed as the convolution (with 
an extra pre-factor of N a /2) of the two following av, l spectral functions for the i = +1 and i = —1 sub-branches, 
respectively, 



B l ^(k',u')= ]T |(0|Fj_ GS , QV , t F p _ h , QV , l Fl GSiQVi jO)| 2 <s(a;'-Z^ l ; QV m Q1/ , t )<S fc , ^, (47) 

J-CPHS-av, l-(C) 

where av = cO, si and Z = ±1. Both the function B 1 q (k',u>') of Eq. (|45|l and the function l|47|l refer to the subspace 
spanned by the tower of states generated by the elementary processes (C) from each initial reduced J-CPHS subspace 
excited state. The summation ^2j^cPHS-au i-(C) on ^ ne right-hand side of Eq. (|4T|I runs over all av — cO, si 
pseudofermion occupancy configurations generated by the small-momentum and low-energy elementary processes (C) 
for the av pseudofermions belonging to the i = sgn(q) 1 = ±1 sub-branch. 

The contributions of the elementary processes (A) to the matrix-element overlaps of the above spectral functions are 
easy to evaluate, and lead to the matrix-element expressions given in Eq. (59) of Ref. [^j . After such contributions are 
accounted for, the problem is reduced to the overlap of the cO (and si) pseudofermion probability amplitudes provided 
in Eq. (69) of that reference. That amplitude is associated with the matrix element (0|Pj_gs, av Fp-h, av PIgs au\ty 
given m Eq. (60) of Ref. [13 and the corresponding av = cO, si pseudofermion spectral function Bq (k',u>') of Eq. 
(|45|l . The generators appearing in the matrix elements given in Eqs. (59) and (60) of Ref. [lOj are expressed in terms 
of the sets of canonical-momentum and bare- momentum values of the distribution function deviations l|13|) - (|17f) in 
Eqs. 1B1(I - I|B5|) of Appendix B. In contrast to the matrix elements of Eq. (59) of that reference, the evaluation of 
the probability amplitude \(0\Fj-Gts,av F p -h, avF_ GS Q „|0)| 2 is a complex problem, which we address in the ensuing 
section. The states Pl GS Q J0) and F^_ h av Fj_ GS Q J0) involved in that probability amplitude describe N° v + 
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AN/// V pseudofermions, whose discrete canonical-momentum values are those of the ground-state and excited-state 
J-CPHS subspace, respectively. Also the state F^_ GS a j0) describes N^ v + AN^/ V pseudofermions, whose discrete 
canonical-momentum values are those of the excited-state reduced J-CPHS subspace. Both the states F_ GS av \0) 
and Fj_ GS al/ \0) refer to densely packed canonical- momentum occupancy configurations. The probability amplitude 

i^ ' = | (0\Fj^cs, av F^ GS al/ \0) | 2 associated with such canonical-momentum densely packed configurations gives the 
av pseudofermion spectral-function lowest-peak weight. By use of the general av — cO, si pseudofermion determinants 
given in Eqs. (71) and (72) of Ref. [ljj, we find the following expression for such a weight, 



4(0,0) 



] \ 2[AT° 1/ +AiV^'„] 



N* 
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x sm 
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x TT TT — 7 

/ = i sm 2 (A^te)AU>(^[^ + 



2JV , 2I 1 " 2 



= cO, si , 



(48) 



where N~®(q) = N® v (q) + AN^(q), N®„(qj) is the ground-state bare- momentum distribution function provided 
in Eqs. (C.1)-(C.3) of Ref. @, and the deviation AN^(q) is given in Eq. Q15|). Moreover, we find that in the 
thermodynamic limit the lowest-peak weight 148fl has the following approximate behavior, 



A< 0,0 > pa TT [l + Ol — )1 • i4<°'°> = 



Here 2A^ is the functional given in Eq. I|4U|I and f av< L reads, 



-1/2+2AJ, 



av = cO, sl ; i = ±1 . 



(49) 



f av ,i = J f{Qau(tqF aJJ ) + sgn(fc)7rj ; f av = JJ t ; az/ = cO, sl , (50) 

where fc is the excited-state momentum relative to the initial ground state, f(Q) — f{—Q) is the function defined in 
Ref. 14j, which appears on the right-hand side of Eq. (24) of that reference, and f av appears in spectral-function 
expressions introduced below. Moreover, S® v is a n, m, and U/t dependent constant such that S° S® 1 — + 1 both for 
U/t — > and for U/t ^ oo and m — ► 0. (From Ref. [lj] we learn that S^g — ► sin(7rn) for U/t — > oo and m — > 0, and 
thus 5gj — > 1/ sin(7rn) in such a limit.) 

The momentum and energy dependence of the spectral-function expression (|44|) is controlled by the convolution 
function, Eq. of the cO and sl pseudofermion spectral functions B l Q a ^(k' ,uj') given in the same equation. The 

main mechanism of such a control is the exotic matrix- element overlap associated with the cO and sl pseudofermion 
elementary processes (C), which generate cO and sl pseudofermion "particles" and "holes" in the vicinity of the Fermi 
points. Although such "particles" and "holes" are not active scattering centers, they are active scatterers whose 
overall phase shifts originate the orthogonality catastrophe leading to the unusual overlaps of the matrix elements 
(0\Fj_GS,av F p -h,av F_ GS av \0). Interestingly, only the overall phase shifts of the scatterers which are not active 
scattering centers and are generated by the elementary processes (C) contribute to the relative weights associated with 
such an orthogonality catastrophe. On the other hand, the value of the corresponding overall phase shifts is solely 
determined by the occupancy configurations of the initial reduced J-CPHS space states generated by the elementary 
processes (A) and (B) from the ground state. Indeed, the latter processes generate the active scattering centers that 
determine the value of the overall phase shifts of the "particles" and "holes" created by the elementary processes (C). 
Therefore, the unusual spectral properties result from the interplay of the elementary processes (A) and (B), which 
generate the active scattering centers, with the elementary processes (C), which generate the active scatterers, as far 
as the overall phase shift associated with the relative weights and the corresponding orthogonality catastrophe that 
controls these properties are concerned. All the av = cO, sl pseudofermions with bare momentum away from the 



15 



Fermi points are also scatterers, but their overall phase shifts only contribute to the lowest-peak weight probability 
amplitude Aa^ given in Eq. (fl%)l . 

Next, let us calculate the convolution function and spectral functions given in Eqs. 1|45J) and (|47(l . The expressions 
of these functions fully define the general spectral functions (|41|l and 144(1 . The canonical-momentum occupancy 

configuration of the state Fj_ GS au \0) associated with the lowest-peak weight probability amplitude is densely 

packed and thus such that {m av ^ + i,m alJ + i) = (0,0). Following the form of the momentum and energy spectra 
provided in Eq. 1)31(1 . the summation on the right-hand side of the Bq^ (fc', u>') expression of Eq. 1(45(1 (and Eq. ((47(0 
over all occupancy configurations generated by the elementary av — cO, si pseudofermion particle-hole processes 
(C) with the same energy and momentum, is equivalent to consider such a summation for all processes with the 
same value of m av . L = [ILj ' Airv al ,](uj' + iv av k') (and m av . L = [IL /2-Kv a v\Lo' = [IlL/2tt] k'). By performing that 

summation we reach the relative weight Aaif* /\(0\Fj^cs, avFp-h,av F~gs aJ0)| 2 corresponding to excited energy 
eigenstates with N^ F L av pseudofermion particle-hole processes. In Appendix B we consider states generated by 
av = cO, si pseudofermion particle-hole processes whose "particle" and "hole" bare-momenta have arbitrary values. 
(The number is not related to the number N^ NF of av = cO, si pseudofermion finite-momentum and finite- 
energy particle- hole processes generated by the elementary processes (A). Instead, it is such that for N a — *• oo its values 
only contribute to the av pseudofermion relative weights when the sum rule = J2i=±i ^m> F i. ls obeyed.) For 
N a — > oo only the relative weights associated with the small- momentum and low-energy av = cO, si pseudofermion 
elementary particle-hole processes (C) in the vicinity of the Fermi points ±qp ai/ lead to finite contributions to the 
av = cO, si pseudofermion spectral function B l Q av (k' , to') of Eq. 1)45(1 . Each of the excited energy eigenstates described 
by the numbers (rn aVt +1, m av , -1) can be multiply degenerate. Therefore, we need to sum up the relative weights 

\N ph ] 

a at f" j Eq. I(B6J) of Appendix B, for all the av pseudofermion particle-hole configurations with the same values of 
(m aVj + \,m aVt Fortunately, for these processes, the general relative weight expressions given in Eqs. I(B6() - ((B9() of 
Appendix B simplify. This follows in part from the symmetry found in Sec. Ill that all the excited energy eigenstates 
generated by the elementary processes (C) from a given initial reduced-subspace state have the same value for the 
four functionals 2A L av defined in Eq. 1(40(1 . We find the following general expression for the relative weights in the 
tower of excited energy eigenstates, 

a a u(m a ^ + i,m a ^ _i) = ^ a av . L (m av , L ) 1 + o( ° j ; av = cO, si . (51) 

L = ±l 

Here a av ^{m a v,i) is the relative weight associated of the spectral function l(4"?jl . It reads, 
a av , b {m av , t ) = @{m av ,,) f[ (2A °" + 3 - 1} = 9(m w , t ) ; i <*v = cO, si ; t = ±l, (52) 

• =1 J 1 \Mctu, t + 1 ) 1 [Zl\ aL ,) 

where T(x) is the usual gamma function. (While, by construction, the integer m a v,c is such that m av . L > 0, for later 
use it is convenient to include the theta-function in the relative weight general expression ((5211 . ) It follows from Eq. 
that, 

a av ,t(l) = 2A^ ; av — cO, si ; 1 = ±1 , (53) 

what reveals that the functional (|12|l is the relative weight associated with the m aUy L = 1 peak of the av, 1 pseud- 
ofermion spectral function 1)47(1 . Moreover, 

Oov(l, 0) = ; Oai,(0, 1) = 2A~l ; av = cO, si . (54) 

The function 1(52(1 has the following asymptotic behavior, 

1 / x2A^-l 

QW, k(m a u, t ) » &(m av , J .,,.., , m av t ; 2A L av ^0; av = cO,sl; t = ±l. (55) 

r (2A^J V J 

Use of the above results in the pseudofermion spectral functions given in Eqs. I(45|) and ((47(1 leads to the following 
expressions for these functions, 
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B l g l av {k',uj r ) = ^2 A av 0) a a y{m a ^ +1 , m aVt _ x ) s(u' - l^- v av ^ m av . L 



m av, + 1, m av, -1 

L 



L = ±l 



l[u>' + iv av k'\ 
4irv ai/ " au .if, """'"V 4nv av /L 



ai?) n <w( 
t=±i 



= cO, si 



(56) 



and 



B l ^(k',Lu') = A av°\ a av, i(m av , L ) s{uj' - l^- v av m a v, S k ,^ 



L 



2nv n 



.4(0.0) a ( )S 



av = cO, ,sl ; i = ±1 , 



(57) 



respectively. In order to distinguish the momentum and energy values of these pscudofcrmion spectral functions from 
those of the spectral function l|41() . k and u, in Eqs. (|45l) . 147fl . (|56() . and 157|l we denote the former momentum and 
energy values by k' and u>', respectively. It follows from Eq. 155fl that for 2A L aL/ ^ the spectral-function expressions 
(|56J) and Ij57(l have the following asymptotic behavior, 



4tt 



f / Zw' \2A' -1 

u eM (^) v*> 



= cO, si ; 
ai/ = cO, si ; t = ±1 , (58) 



for small finite values of 1[lu' + iv av k'] and Ilu' , respectively. 

Let us consider the convolution function (|45|l for vanishing energy values of the order of 1/L. For instance, the 
six smallest discrete energy values read IAlu = 0, IAlu = [2tt/L]v s i, IAlu = [2ir/L]v c o, IAlu = [4tt/L]v s i, IAlo = 
[2tt/L][vco + v s i], and IAlu = [Att/L]v c q such that IAlu < [An / L]v c q. For vanishing energy values we define that 
convolution function as B 1 ' 1 (Alu, v) = B 1 ' 1 (Alu, Ak), where by use of the first expression of Eq. I|45|l we find, 



L? M (Aw,Afc) = 2tt 



N„ 



Co 



D S(ALu)S(Ak) ; D Q 



n 



au— cO, s 



\2A„„ ' 



IAlu = ; ZAfc = : 



Co 



B l ' l {Au>, Ak) = 2ir [ — ) D a 5(lAu - u(C)) 5(lAk - fc(C)) ; IAlu = lu(C) ; v = 



fe(C) 



(59) 



Here Co an d 2A Q „ are the functionals given in Eq. (|39[1 and for the above five energy values such that IAlu > 
the relative weight a reads: ao — 2A L av for lu(C) = \p,Tr/L]v av , k(C) = i[2n/L], and v — to al/ ; a — 2A L av 2A~„ 
for w(C) - [47r/L]« atf , fc(C) - 0, and l/v = 0; a = 2A^ V 2A^ C for = [27r/L][t; c o + v s i], fc(C) = t[47r/L], 

and = l [v c o + v s i}/2 where we introduced the index av such that cO = si and si = cO; ao = 2A^„ 2Agp for 
lu(C) = [2tt/L][v c0 + v s1 ], k(C) = 0, and l/v = 0; and a = 2A^[2A^ + l]/2 for uu(C) = [47r/L]v av , k{C) = l[4tt/L}, 
and v = Lv au . All these expressions refer to av = cO, si and l = ±1. Similar ao expressions can be derived for larger 
energy values of order 1/L. 

In turn, for small finite values of IAlu we use the av pseudofermion spectral function expression (|56|l and the 
variables x = oj'/Alu and y = sgn(u)fc' [v c q/Alu] where Alu = (10 — IAE) in the second expression of Eq. l(4"5|l and find, 
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2ttv c0 \4ir y 'v c0 v s /LJ J 
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(60) 



The expressions given in Eqs. (|59|l and (|6U|I combined with Eqs. I|41|l and 144|l provide the desired general spectral 
function expression, except for k and lu values such that lu w LV av (k — Ika) + Iluq, where av — cO, si and t = ±1. 
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In the ensuing section we use the function (|44|l to derive closed-form analytical expressions for the finite-energy one- 
and two-electron spectral-weight distributions. 

For small finite values of IAlu = l(ui — IAE), use of the relative weight asymptotic expression (|55|l in Eq. (|60|l for 
the whole x and y integration domains leads to the following asymptotic behavior for the convolution function (|60|l . 



m^v)~/^e(lA.)(-±p=y 2+C0 ; 1 = 0,1,2,..., l = ±l. (61) 
Here Co is the functional given in Eq. (|39[) and the function F (z) reads, 
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Vsl 
V C 



r(2A^) 

2At, -1 



r(2A^ ) 

x + sgn(z) l y 



(62) 



where 2A s i is the functional given in Eq. (|39(1 and Do is defined in Eq. 1)59(1 . Note that both the values of the functional 
Co = Co({qi},{q'i},{qj},Wj}) defmed in E q- (EH and of the function F (z) = C {z, {<?,:}, {q'J, {q^}, {</,}) of Eq. 
(16211 depend on the set of bare-momentum values {qt}, {q'i}, {qj}, W j} which define the active-scattering-centers occu- 
pancy configurations of the initial reduced-subspace excited energy eigenstates generated by the elementary processes 
(A). Such a dependence follows from the functional character of the quantity 2A^„ = 2A L a „({qi}, {q'i}, {qj}, Wj}) 
defined by Eq. (|40|l . The dependence of the value of the latter functional on the set of bare-momentum values 

Ui}, Wi}, {qj}, {q'j} occurs through the scattering phase shift Q%i N ■ \iq Fav , {qi}, {q' t }, {q 3 }, {q'j})/2 appearing in 
Eq. (|40|) . which is defined in Eq. (|36() . Thus, we should keep in mind that the values of the functionals 2A^„ and Co 
appearing in the spectral-function expressions derived below are different for each pseudofcrmion bare-momentum oc- 
cupancy configuration defining a specific excited energy eigenstate generated from the ground state by the elementary 
processes (A) and (B). 



V. WEIGHT DISTRIBUTION EXPRESSIONS GENERATED BY THE DOMINANT PROCESSES 



In this section we derive analytical closed-form expressions for the one- and two-electron weight distributions 
associated with the function (|4"Tll for (k, w)-plane regions whose weight is generated by the dominant processes. 
These processes are characterized below in terms of rotated-electron processes. Our study includes the derivation of 
the weight distribution expressions valid for k and to values such that u> ~ i v av {k — Iko) + Iloq, where av = cO, si and 
i = ±1. Some of the analytical expressions found below refer to the vicinity of the weight distribution singularities 
and edges. Importantly, the spectral features observed in real experiments correspond to such singularities [3, 0? ] • 



A. THE DOMINANT PROCESSES AND THE FUNCTIONAL 3 {k, ui) 



The spectral weight distribution associated with the terms of the general finite-energy spectral function (|41l) of 
index i > decreases very rapidly for increasing values of the number i of extra rotated-electron pairs. Indeed, for 
the excited energy eigenstates associated with the spectral-function contributions of increasing i value the functional 
Co of Eq. (|39|) also has increasingly larger values. Typically, the contributions of orders i larger than i — 1 are 
beyond numerical measurability, once the spectral weight associated with the i = term of the spectral function 
(|4H corresponds ingeneral to over 99% of the whole spectral weight. (This is confirmed for the one-electron spectral 
function in Ref. [jjj-) Thus, for practical applications it is in general enough to consider the elementary processes 
associated with the i = term only. 

The smallness of the weight associated with the i > terms of the spectral function (|41|l results from the form 
of the function l|44|) . Indeed, for increasing values of the index i the expression of the operator Q 1 ^ ' av associated 
with the operator B^^' 1 , Eq. (58) of Ref. 01 , has in terms of rotated-electron operators an increasing number of 
rotated-electron particle-hole pair operators |10| . Thus, for increasing values of i application of such an operator onto 
the ground state produces elementary processes (A) and (B) which generate an increasing number of active scattering 
centers. Since the value of the functional Co given in Eq. (|39f) is in general an increasing function of the number of 
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generated active scattering centers, it follows that the i > contributions to the spectral function are very small. 
Below it is confirmed that most singularities of the latter function occur at (fc, cj)-plane isolated points and lines, are 
of power-law form, and are controlled by processes which generate excited states whose Co values are smaller than 
two and one, respectively. Thus, a quite good approximation corresponds to replacing the general spectral-function 
expression (|41H by its i — term. The excited-energy-eigenstate deviations generated by the dominant processes 
associated with the i — operator 0jy o ■ on the right-hand side of Eq. (32) of Ref. 0] obey the two selection 
rules given in Eq. (41) of the same reference. The numbers l — A/j? ; in that selection rule are those of Eq. 

(30) of Ref. 0] specific to the corresponding operator 8^- ■ defined in Eqs. (27) and (28) of the same reference. 
Furthermore, all excited energy eigenstates whose deviations do not obey that selection rule are generated by the 
i > operators on the right-hand side of the latter equation. It follows that a good approximation for the spectral 
function consists in using expression 141fl for i = and the summation on the right-hand side of Ecj_. (|41|> over the 
J-CPHS subspaces whose number deviations obey the sum rules (18), (19), (39), and (43) of Ref. [Hj and selection 
rules given in Eqs. (21) and (41) of the same reference. 

Moreover, for the one- and two-electron spectral functions nearly the whole weight is concentrated in the (fc, w)- 
plane regions associated with creation by the elementary processes (A) of none, one, and two active scattering centers 
away from the av — cO, si Fermi points and av ^ cO, si limiting bare-momentum values. Concerning singular 
weight features, we find below that for finite values of the on-site repulsion U only exceptionally and for some spectral 
functions and isolated (k, w)-plane points are these features of (5-function type. In general, such singularities are 
instead of power-law shape. Power-law spectral-function behavior occurs for (k, w)-plane regions in the vicinity of 
points or lines corresponding to different reduced J-CPHS subspaces than that associated with such regions. In order 
to describe such an effect it is convenient to associate each (fc, w)-plane point with the minimum value of the functional 
Co, Eq. (32) ; of the excited energy eigenstates generated by the elementary processes (A) and (B) that correspond to 
such a point. We thus introduce the functional, 

3 (fc, w) = minCo • (63) 

In general, for the same regions of the (fc, w)-plane the spectral function (|41|) has contributions from different func- 
tions B l ' l (k, us), associated with different J-CPHS subspaces of the same CPHS ensemble subspace. The value of the 
functional (|63[) corresponds to the smallest value of Co of all states associated with the (k, w)-plane point into consid- 
eration. For two-dimensional (k, w)-plane regions associated with a given reduced J-CPHS subspace, the value of the 
functional £o> Eq. I|39|l . is in general a continuous and smooth function of k and u). Therefore, in that case the value 
of the functional 3o(fc, oj) is also in general a continuous and smooth function of k and u>. Importantly, the spectral 
weight power-law features correspond to the points or lines where ^o(k, ui) shows discontinuities. (The inverse is not 
always true.) Typically such discontinuities occur in the border lines of the (fc, <x>)-plane domains associated with the 
reduced J-CPHS subspaces. Such subspaces can correspond to a {k, w)-plane isolated point, line, or two-dimensional 
domain. In general, the corresponding value of the the functional Co, Eq. (|39|l . is smallest for the point-subspaces 
whose deviation numbers and numbers are given in Eq. (^J. It is also in general smaller for the line-like reduced 
J-CPHS subspaces than for the subspaces corresponding to a (k, w)-plane two-dimensional domain. Therefore, when 
the reduced J-CPHS subspace (k, cj)-plane domains correspond to isolated points or lines, the discontinuities of the 
value of the functional (|63|) occur at these points or lines. It also shows discontinuities at the (k, cj)-plane border lines 
defining the limits of the reduced J-CPHS subspace two-dimensional domains. 

The reason why the spectral function has power-law behavior in the vicinity of the (k, w)-plane isolated points and 
lines, where the value of the functional (|63|l has discontinuities, is the small energy gap between these regions and the 
point or line corresponding to the discontinuities. The value of the functional (|39|l is in general different at the point 
or line than in its vicinity. Moreover, in general that value is smaller for the former line or point than in its proximity. 
Thus, the small gap gives the minimal excitation energy of the tower states associated with the (fc, w)-plane point 
where one is calculating the spectral function and which are generated by the elementary processes (C) relative to that 
of the initial states corresponding to the point or line where the value of the functional l|t)3|) has a discontinuity. We 
recall that all states generated by the elementary processes (C) have the same value for the functional Co, Eq. (|39|) . as 
the initial state, which corresponds to the point or line where the value of the functional i|63|) has a discontinuity. The 
occurrence of the small gap implies that the corresponding convolution function l|60|) has the power-law asymptotic 
form 1|61[) . Note that such a power-law behavior is controlled by the exponent — 2 + Co, which is negative provided that 
Co < 2. Furthermore, once the value of the functional (|39|l is smaller for these excited states than for states generated 
by elementary processes (C) from initial states corresponding to (fc, o;)-plane points located closer to the point where 
we are calculating the spectral function, the spectral function itself has a power-law shape, as confirmed below. 

In turn, if in the vicinity of a (fc, w)-plane point the value of the functional I|tj3f) has no discontinuities, the spectral 
function expression l|44|l at that point involves integrations over gapless contributions from transitions generated by 
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the elementary processes (C) from initial states corresponding to a small two-dimensional (k, w)-plane domain just 
below (I = +1) or above (I = —1) the point. Once the value of the functional Ho(k, u>) is in this case a continuous 
and smooth function of k and u>, its value remains the same in the infinitesimal vicinity of the point. In this case the 
spectral function has no power-law behavior, as confirmed below. The lack of such a behavior follows from the gapless 
character of the excitations contributing to the spectral function at that point. Below we consider the spectral- weight 
distribution behavior generated by creation of none, one, and two active scattering centers away from the av = cO, si 
Fermi points and av ^ cO, si limiting bare-momentum values. For one- and two-electron spectral functions the weight 
generated by processes involving the creation of more than two active pseudofermion or hole scattering centers away 
from the av — cO, si Fermi points and av ^ cO, si limiting bare-momentum values is very small and can in general 
be neglected. 

Our results reveal that the functional (|63|1 plays a major role in the control of the spectral weight distribution. 
Increasing the number of created active pseudofermion and pseudofermion hole scattering centers increases in general 
the value of that functional. Moreover, the summations in expression l|44|) over the states generated by the elementary 
processes (A) further increase the value of the power-law exponent — 2 + Co appearing in the convolution function 
B61J1 . Indeed, these summations add integer numbers to the exponent —2 + (q: creation of each active scattering 
center adds 1 to that exponent, as confirmed below. A central point is that most spectral weight singularities are 
of power-law type and correspond to the discontinuities of the value of the functional (|63|l in the (k, w)-plane. Also 
the weight associated with creation of two or more active scattering centers away from the av = cO, si Fermi points 
and av ^ cO, si limiting bare-momentum values is controlled by that functional. Furthermore, the same functional 
controls the exceptional occurrence of <5-function singularities: such singularities corresponds to the (k, w)-plane points 
where its value vanishes. 

Each spectral-weight distribution (|44ll corresponds to a J-CPHS subspace. We start by considering the spectral- 
weight distributions generated by creation of two active pseudofermions and/or pseudofermion holes scattering centers 
away from the av — cO, si Fermi points and av ^ cO, si limiting bare-momentum values. Thereafter, we consider 
spectral-weight distributions (|44|l associated with contributions from processes involving creation of none and one 
active scattering center away from the av = cO, si Fermi points and av ^ cO, si limiting bare-momentum values. We 
find below that the latter processes lead to important point- and line-like singular spectral-weight features. 



B. TWO-DIMENSIONAL (k, w)-PLANE REGIONS WHERE THE SPECTRAL WEIGHT IS 
CONTROLLED BY TWO ACTIVE-SCATTERING-CENTER CREATION 

The parametric equations that define the two-dimensional (k, w)-plane domains generated by creation of two active 
scattering centers away from the av = cO, si Fermi points and av ^ cO, si limiting bare-momentum values are of the 
general form, 

k = lk a „ t (q, q')\ uj = luj alJ: a / v > (q, q') . (64) 

For instance, 

k = l[ko + ci q + c[ q'] ; w = lu> av , a 'u'(q, q') = l[^o + ci e au (q) + c[ e a / v >(q')] , av = cO, si , aV = cO, si ; 

k = l[k + ci q + k cv - q] ; uj = luj av<a > v >(q, q') = l[u> + c\ e au (q) + el v ,(q')] , av = cO, si , v' > ; 

k = l[k + ci q + q'] ; w = luJ av<a ' v '(q, q') = l[wo + ci e au (q) + e° su ,(q')] , av = cO, si , v' > 1 . (65) 

Here the index / — ±1 is that of the corresponding spectral function B l j^(k, u>) of Eq. Q4ip. the momentum k CV) energy 
ujq, and momentum fco are given in Eqs. (|27|) . i|32|) . an d 134|) . respectively, and the constant c\ [and c'J is such that 
ci = +1 (and ci = —1) for creation of a av pseudofermion (and a av pseudofermion hole). For simplicity, in Eq. 
(|65|l we considered that the two created active scattering centers are both av = cO, si pseudofermions and/or holes, 
a av = cO, si pseudofermion or hole and a cv ^ cO pseudofermion, and a av = cO, si pseudofermion or hole and a 
sv 7^ si pseudofermion. Generalization to other active-scattering-center choices is straightforward. 

The summation over the excited energy eigenstates generated by the elementary processes (A) which contribute to 
the spectral weight at the (k, w)-plane point defined by Eq. I|65l) is performed in Appendix B. The momentum and 
energy of these states relative to the ground state is close to k and u, respectively. These contributions occur through 
generation from these states, by the elementary processes (C), of new excited energy eigenstates whose momentum 
and energy relative to the initial ground state are precisely k and uj, respectively. Since the former excited energy 
eigenstates generated by the elementary processes (A) have very close momentum and energy, except for vanishing 
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contributions of order 1/N a the corresponding value of the functional H63|) is the same. For finite values of that 
functional it is found in Appendix B that the spectral function has following general form, 



7rC c C s 



+l/v st 

dz Fq(z) 

■l/v.i 



Q, \ Co (?,?') ^v cQ v sl 



v s i I Co(q, q') \v av {q) - v a > v >{q')\ 



l = ±l. (66) 



While this expression and all other spectral-weight distribution expressions obtained in this paper correspond to 
densities < n < 1 and < to < n, taking the limit to — > in such expressions leads to the correct m = results 
[T^ |. The expression Ij66(l is not valid in the vicinity of isolated (k, w)-plane points whose weight is associated with 
point-subspaces and in the proximity of the branch lines studied below. In the vicinity of these points and lines one 
must consider the contributions from processes involving the creation of none and one active scattering center away 
from the av = cO, si Fermi points and av ^ cO, si limiting bare-momentum values, respectively. On the right-hand 
side of Eq. (|66|l . Co = Co(<2S q') is the value of the general functional (|39Jl specific to the two active-scattering-center 
excitation under consideration, F$(z) is the function 1)62(1 . and the small but finite energy fl corresponds to the range 
of the elementary processes (C). We recall that fi is treated as a parameter whose value for fixed finite values of n, to, 
and U/t is determined by imposing the exact sum rule to the full spectral-function expression in the end of the 
calculations. Q also appears in the spectral- function contributions from none and one active scattering center creation 
away from the av = cO, si Fermi points and av ^ cO, si limiting bare-momentum values, as confirmed below. The 
general spectral- function expression (|66|l is valid for finite values of the functional l|63|) . For regions of the (fc, w)-plane 
where the value of that functional vanishes the energy f2 is such that 17 — > and the spectral-function expression has 
a different form, as discussed below. 

The k and u> dependence of the weight-distribution expression ((66JI occurs mainly through the q and q' dependence 
of the value of the functional Co(os '/) and through the absolute value of v a „(q) — v a i u i(q'). The created active- 
scattering-center bare momentum values q and q' are related to the weight-distribution value of k and lo by Eq. I|65[) . 
We emphasize that there are two types of (fc, w)-plane points: (i) those associated with one pair of (q, q') values and 
(ii) those corresponding to two different pairs of such bare-momentum values. In the case (ii) the spectral function 
at the (k, w)-plane point is given by the sum of two terms of the general form given in Eq. (|66|l . The shape of 
the two-dimensional sub-domains of type (i) and (ii) depends on the specific form of the parametric equations iJBSJ. 
While the goal of this paper is the derivation of general analytical expressions for the one- and two-electron spectral- 
weight distributions of the metallic phase, examples of such sub-domains will be given elsewhere for specific one- and 
two-electron spectral functions. 

The two-active-scattering-center (k, w)-plane two-dimensional domain is bounded by well defined lines. Also the 
two sub-domains (i) and (ii) are bounded by well defined lines. An important general property is that the part of 
the limiting line of the two-dimensional domain which is also a limiting line for the sub-domain (ii) [this excludes the 
limiting line between the sub-domains (i) and (ii)] is such that v av (q) — v a ' u '(q'). We call such a line border line. The 
shape u> = u>BL,(k) of a border line is defined by the following parametric equations, 



^Bh{k) = 1[lu + ci e av (q) + c[ e a >„>(q')] 8v a »(q),v a ,„,(q>) 5 k = lk au , a ' v >(q, q') $ Vau (q), v a , u ,(q') ■ (67) 

According to expression (|66|l . the weight distribution has singular behavior at that line. Note that such an expression 
defines the k and u> dependence of the weight distribution for the regions just below (I = +1) or above (/ = —1) 
the border line. We recall that the border line corresponds to a discontinuity in the values of the functional 163fl . 
Thus, the weight distribution has power-law behavior in the region just above (I — +1) or below (I = —1) the border 
line. The k and ui dependence of the weight distribution in that region is not given by Eq. I|66() . The study of 
that dependence proceeds much as for the branch lines considered below and will be carried out for specific weight 
distributions elsewhere. The same applies to the region just above (/ = +1) or below (/ = —1) the remaining limiting 
lines of the two-dimensional two-active-scattering-center (k, w)-plane domain. 

We emphasize that depending on the specific weight distribution, the dominant processes can correspond to different 
choices of active-scattering-center pairs. Thus, it can occur that the same (fc, cj)-plane region involves contributions 
from several expressions of general form given in Eq. (|66|) . In that case the weight distribution at a (k, w)-plane point 
is given by the sum of several terms of the general form provided in that equation. Moreover, one must also add the 
spectral features generated by creation of none and one active scattering center away from the av — cO, si Fermi 
points and av ^ cO, si limiting bare-momentum values, which we derive below. 
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C. THE POINT-LIKE AND LINE-LIKE SPECTRAL- WEIGHT DISTRIBUTION FEATURES 



The smallest reduced subspaces are point-subspaces whose deviation numbers and numbers are given in Eq. 
Analysis of the general spectral-function expression (|41|) for i = reveals that the most divergent singular power- 
law spectral features appear in the vicinity of the (fc, w)-plane points, (lk Q , Iloq), corresponding to such subspaces. 
Such a weight results from processes that do not involve creation away from the av — cO, si Fermi points and 
av 7^ cO, si limiting bare-momentum values of active scattering centers. For (fc, w)-plane regions in the vicinity of the 
corresponding (Iko, Iloq) points the dominant contributions to the general expression l|44J) are generated by elementary 
processes (C) whose initial state corresponds to the point-subspace. Here ujq and fco are given in Eqs. and J2IJ), 
respectively. Such contributions lead to B l h f(k, uj) « 0(£1 — l\u> — Iuiq}) B 1,0 (uj — Iuiq, [uj ~ luio/k — Iko]) and thus the 
weight distribution has the following power-law behavior, 



where I = ±1 and the velocity v = [uj — loJo]/[k — Iko] is always such that \v\ > v s % and v ^ ±v c o- Depending on 
the specific point-subspace and weight distribution, the domain of z = 1/v values corresponding to the small finite 
spectral-weight (fc, w)-plane region where the expression (|68|) is valid is in general bounded by two of the four values 
— 1/ U c0> — and 1/v c q. This expression refers to small finite values of [uj — luJo], the functional £o and the 

function F (z) are given in Eqs. and respectively, and the value of the corresponding functional 2A^ ) 1 and 
2Afi whose general expression is given in Eq. (|40|l is that of the initial point-subspace. It is not valid for k and u> 
values such that uj ~ t v av (k — Iko) + Iujq, where av — cO, si, l = ±1. (The expression for such k and uj values is given 
below.) When the exponent — 2 + £ is negative, the weight distribution IpjSI) has a singularity at the (lk , Iuj ) point. 
In general, the value of that exponent increases for increasing value of |fco|. Thus, it is larger for the spectral function 
expression in the vicinity of the (Iko, Iujq) points corresponding to the second point-subspaces than in the vicinity of 
the points associated with the first point-subspaces. 

For the small (fc, w)-plane region in the vicinity of the point (Iko, ^o)j where the expression (|68|) is valid, the 
contributions from functions of form (|44|) associated with J-CPHS subspaces other than the point-subspace are very 
small. Therefore, in that region the spectral function l|41(l has for i = approximately the form given in Eq. (|68|l . 

Let us now consider the contributions from processes involving the creation of a single active scattering center away 
from the av — cO, si Fermi points and av ^ cO, si limiting bare-momentum values. These contributions lead to the 
next more divergent singular features, which arise in the vicinity of branch lines. In this case the energy eigenstates 
which span a reduced J-CPHS subspace are associated with a (fc, w)-plane line that corresponds to a discontinuity in 
the values of the functional l|63f) . By changing the bare momentum q of the active av pseudofermion or pseudofermion 
hole scattering center, one generates a branch line in the (fc, w)-plane. Often one of the end points (or both end points) 
of such a line coincides with a (Iko, l^o) point associated with a point-subspace. For each CPHS ensemble subspace, 
at least one of the end points of the most divergent of such branch lines coincides with a first point-subspace point 
(Iko, 1<jJq)- The active av pseudofermion (and av pseudofermion hole) scattering center can belong to any branch 
(to the av = cO, si branches). A av — cO, si branch line is generated by creation of the av pseudofermion (or av 
pseudofermion hole) for all available bare-momentum values in the domain \q\ s [qp a ^, q% v ] (or |g| £ [0, g^,,]). 

The parametric equations that define the (fc, w)-plane points belonging to a av pseudofermion (or av pseudofermion 
hole for the av = cO, si bands) branch line is of the general form, 



fc = lk av (q) = l[ko + ci q] ; uj = luj av (q) = I [uj + Ci e av (q)] , av = c0,sl; 

k = lk cv (q) = l[k + k cv - q]\ uj = luJ cv (q) = 1[ujq + e%(q)] , cv^cO; 

k = lk sv (q) = l[k + q] ; uj = luj su (q) = 1[uj + e° sl/ (q)] , sv^sl. (69) 

Here the index / = ±1 is that of the corresponding spectral function B l j^(k, uj) of Eq. 141fl . the momentum fc c „, energy 
ujo, and momentum fco are given in Eqs. I|27l) . I|32l) . and (|34fl . respectively, and the constant c\ is such that c\ = +1 
(and c\ = —1) for creation of a av pseudofermion (and a av pseudofermion hole). 

Let us consider a (fc, w)-plane point located just above (I = +1) or below (/ = — 1) the branch line whose momentum 
fc obeys the relation given in Eq. I|tj9fl and the energy uj is such that l(u — luj aiy (q)) is small and positive. The weight- 
distribution expression at that point is controlled by the elementary processes (C), which generate from the initial 
excited energy eigenstates corresponding to the line a set of tower states whose momentum and energy relative to 
the ground state are precisely fc and uj. By performing the summations of the general expression (|44|l over the initial 
excited energy eigenstates generated from the ground state by the elementary processes (A) which correspond to line 
points in the vicinity of such a point, in Appendix B we derive the following spectral-weight distribution expression, 
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B l „(k, w) « e(n - l[u - lco av (q)}) Q(l[w - ^av(g)]) C au (q) ^J^^ ] ) ; Co(<z)>0, (70) 



where 



sgn(g) c aw ; / Z[w - Zw a „(g<)] \ 



=sgn(q) c aL , 



This expression does not apply to k and lo values such that to ~ LV av (k—lko)+lu>o, where av = cO, si and i — ±1. The 
functional Co = Co(?) appearing in expression l|71|) is given in Eqs. 1391) and Fq(z) is the function 1)6211 . Furthermore, 
Can — 1 for av — cO, s^, c Q „ = — 1 for <w = cv ^ cO. 

The unity added to the exponent — 2 + Co of the convolution function 1)61)1 and spectral function expression 1)680 to 
reach the branch-line exponent — 1 + CoGz) 01 expression (|70|) . stems from the summations performed in Appendix B 
over the suitable initial excited energy eigenstates generated by the elementary processes (A). For a singular av branch 
line the corresponding power-law exponent — 1 + Co (?) is negative. Since when — 1 + Co (<?) < the weight distribution 
shows singularities at the branch line, we expect that in this case spectral peaks will be observed in experiments. 
This is confirmed for one-electron removal in Ref. [? ] for both the TTF and TCNQ spectral dispersions of the 
quasi-lD organic compound TTF-TCNQ, following the preliminary results of Refs. 0, f° r the TCNQ stacks of 
molecules. When the region above (/ = +1) or below (I = —1) the branch line is not contained in a (fc, £x/)-plane 
domain associated with a two-active-scattering-center reduced J-CPHS subspace, the expression (|70|l describes the 
spectral- function behavior when — 1 + Co(<?) < 0> whereas < — 1 + Co(<?) < 1 is a sign of near absence of spectral 
weight. In contrast, if that region is contained in such a domain, the expression (|70() describes a spectral-function 
weight edge when < — 1 + Co(<z) < 1- Finally, — 1 + Co (9) > 1 is always a sign of near absence of spectral weight in 
the vicinity of the line. 

For momentum values k and energies lo such that u> ~ t v au (k—lko)+ltjjQ, where av = cO, si and l — ±1, the spectral 
function corresponds to the vicinity of a av = cO, si branch line end point. However, in this case the expression 
of the spectral function in the vicinity of the branch-line is not that of Eq. 1)70)1 . This results from a resonance 
effect: the branch line group velocity v av (q) equals the velocity iv av associated with the av, 1 pseudofermion particle- 
hole excitation sub-branch generated by the elementary processes (C). Moreover, the av = cO, si branch line is in 
this limiting situation generated by these elementary processes. Indeed, one can consider that in a first step the 
branch-line av — cO, si pseudofermion or hole is created at the q — iq® Fav Fermi point by the elementary processes 
(B). This first step corresponds to a given reduced J-CPHS subspace. One then accesses the corresponding J-CPHS 
subspace by moving the av — cO, si pseudofermion or hole along its energy and momentum dispersion. This second 
small-momentum and low-energy process does not involve creation of any quantum object and thus corresponds to an 
elementary process (C). When moving the av = cO, si pseudofermion or hole in the vicinity of the q = Lq%* au Fermi 
point, its energy dispersion is linear in the momentum. It follows that for these values of k and u> the elementary 
processes (C) replace the elementary processes (A) in what the generation of the av branch line is concerned. While 
the elementary av = cO, si pseudofermion particle-hole processes (C) generate the branch line, the elementary av 
pseudofermion particle-hole processes (C) play the usual role of these processes in the general convolution function 
(145)1 . (As in Sec. IV, here we have used the index av such that cO = si and si = cO.) It follows that for values of 
k and lo obeying the relation lo iv av {k — Iko) + Iloq, where 1 — ±1, the convolution function associated with the 
av = cO, si spectral function is given by, 



1 pAuj pLljj'/v al , 

B l f„ t (Aw,u) « — / du' dk'B l A° (ALo/v-k',ALo-uj')B l A (k^Lo'); I = ±1 , (72) 

2tT Jo J llo> /v au -iAqu n * V ' V ' 

where Aqu ne = l[2ir/L]x such that < x < 1 is the effective momentum width of the branch line. Here we 
consider the more general case of the functional 2A~£ having values in the domain < 2A Q ;£ < 1 and define x 
as x — [l/N a \( 1 / 2A °'v~ 1 \ The introduction of the branch-line momentum effective width Aqu ne provides a good 
approximation for the value for the power-law multiplicative constant given below in Eq. 1)74)1 . In turn, the power-law 
k and lo dependence obtained below is exact. The function (|72)l differs from the general convolution function l|43|) 
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in the momentum integration domain. Indeed, now the role of one of the elementary av = cO, si pseudofermion 
particle-hole processes (C) is to generate the av branch line. 

By use of the same methods as for the general spectral function expression (|70|l in the vicinity of a branch line, in 
Appendix B it is found that for k and u> values such that u> ~ iv au (k — Ika) + Iluq the spectral- weight distribution 
expression is given by, 

B l M (k,uj) « - l(w - + LV al/ (k- lk )])\ Q[l(u - l[u) + tv av (k - Zfco)])J 
v n ( Ku + i v av (fc - lk )]) \ -!+cL, . 

x G av A — ; = cO, si ; t = ±1 ; I = ±1 , (73) 

V 47Ty^cO v sl J 

for £° t > 0, where the exponent £° t is given in Eq. HB20(I of Appendix B, 



il < / l(u> - l[oj + LV av (k - lk )]) 



Cctv, l — ~ — ■p^—pr\Q\V av — V & v\^ 
Z7T U r O, I. V 





+ 6[u a5 [l \-v av ) l dz\- ■ — — — , (74) 

and the function F® v L (z) is provided in Eq. 1|B21(1 of that Appendix. This weight-distribution expression refers to the 
proximity of a av = cO, si pseudofermion or hole branch-line end point and differs from the corresponding general 
spectral function expression given in Eq. (|70|l . The latter expression applies in the vicinity of the av — cO, si branch 
lines when the branch- line group velocity v au (q) is such that v av (q) ^ i v av for i = ±1. As v av (q) — > c v au , the spectral 
function k and u> values correspond to the vicinity of a branch-line end point and thus it is given by expression l)73|) 
instead of 1)70(1 . Moreover, in the proximity of the (Iko, lu>o) point but for values of k and u> corresponding to the 
region where the expression (|68[) is valid, the spectral function is given by that expression. There are very small 
cross-over regions between the (k, w)-plane regions where these different expressions apply. 

D. THE <5-FUNCTION SINGULARITIES 

For some spectral- weight distributions there are isolated (fc, w)-plane points, (Iko, luio), associated with point- 
subspaces such that AN au — for av = cO, si, N av = for av ^ cO, si, and L c _ 1 / 2 > and/or L s ^ _ 1 / 2 > 0. These 
point-subspaces correspond to (Iko, lu>o) points such that, 



l(»-l[u[) + .» m (l=-lt )l) 

n 



dz 



k = 7T i c , -1/2 I ^0 



2 M^c, -1/2 + 2(J,qH L s< _ 1 / 2 



L c.-i/2 and/or £ s ,_i/ 2 >0, (75) 



where the energies 2/x and 2^lqH are given in Eq. I|21|) . Interestingly, the four functionals 2A L au of Eq. I|4(J|) and thus 
the value of the functional (|63|l vanish at these points for all values of U /t. In this case the expression of the general 
operator 0^ k defined in Eq. (27) of Ref. [T(i| does not include any pseudofermion operators, since no pseudofermions 
are created nor annihilated under the ground-state - excited-energy transition. Moreover, the elementary processes 
(C) give no contribution, since when the four functionals 2A L a[/ of Eq. I|4UI) vanish the energy which defines the 
energy range of these processes also vanishes. In this case the operator QKj- k is such that in the matrix element 
(/; CPHS.L\Q l N k \GS) of Eq. (26) of Ref. [l^ the energy eigenstate \f; CPHS.L) is the ground state itself. Thus, 
for such point subspaces the matrix elements derived in Sec. V of Ref. must be replaced by (GS\<d L j^ k \GS) and 
the spectral function has at the corresponding (fc, w)-plane points the following 5-function behavior, 



\{GS\@{f h \GS)\ 2 / \ 
Bjv(fc,w) = - ~cU' — 6 {^- l ^)S k , lko ; l = ±l. (76) 



The (5-peak singularities of Eq. I)76|) correspond to isolated (Iko, l^o) points and occur for all finite values of U/t. 
Another example of <5-peak singularities occurs for the one-electron spectral- weight distributions. Such a behavior 
occurs for some cO or si branch lines (JBS}, in whose vicinity the spectral function has the form (|TU1) for U/t > 0, and 
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for U/t — > the value of the functional (|tj3|) behaves as 3o(k, u>) = 3o(fc, u> ai/ (q)) — + 0. Also the overall phase shifts 
Qauii'Qpav)/'^ behave as Q av {iq Fav ) /2 — ► for = cO or = si and i = ±1. In this case the spectral-function 
expressions derived in this paper remain valid, provided that one considers in these expressions the limit 2A^„ — > 
for the four functionals given in Eq. I|40(l . The energy f2 which defines the energy range of the elementary processes 
(C) behaves in this case as ft — > for ~2o(k, oj) — > and, therefore, the weight distribution becomes 5-function like. 
The function l(5U)l is such that f aV)l — \J f{— n) — > 1 as Qaui^Fav) — * 0. Indeed, /(Q) = f(—Q) is the function on 
the right-hand side of Eq. (24) of Ref. 14] such that f(n) = f(-n) = 1. Moreover, 5^5^ -> 1 as U/t -> and 
thus the <5-peak weig ht 27r(l/JV ) Co D of Eq. © is such that (l/N a )^D -> 27r as U/t — ► 0, once for the above 
branch lines also Co ~^ 0. It follows that in this case the convolution function (|45|l has not the form given in Eq. Q6ip. 
Instead, as [//< — > it is given by the first expression of Eq. (|59jl and the asymptotic expression H70j) is not valid for 
U/t — > 0. By performing calculations similar to those of Appendix B but with the convolution function given by the 
first expression of Eq. I|59|) , instead of its asymptotic expression ()61|l , one arrives to the following expression for the 
spectral-weight distribution at the branch line, 



= ^ / dq' (Jj-^) ° S(w-lu, au (q'))s(k-lk av (q')) 

= ^-5{u-lu av {q)); Co(g') ^ 0; l = ±l, (77) 

where the functions k av {q) and oj au (q) are defined in Eq. i|(j9fl an d the q' integration is over the branch line points 
such that CoW) — ► 0. 

For two-electron spectral-weight distributions a similar behavior occurs for the weight generated by two active- 
scattering-center creation away from the otv = cO, si Fermi points and av ^ cO, si limiting bare-momentum values in 
two-dimensional regions of the (fc, w)-plane where 3o(fc, oj) — > 0. Again, for such regions the convolution function 145|) 
has not the form given in Eq. (|61|l but that corresponding to the first expression of Eq. I|59|l . Thus, as Ho(k, oj) — > 
the expression l|66|l is not valid and instead the spectral-weight distribution has the following behavior, 



D f f ( 1 \^°( 9 ' 9 ') / \ / \ 

B hi k iu) = 2 ^ c ° c dq dq' (^f) S\u - lu av , a > v >{q, q')j S[k - lk av , a > v >(q, q')j 

dq'5(u)-lu av , alv ,{q{k,q'),q')y, Co(?. ?')^0; I = ±1 . (78) 



2n C c C s 

i(q, q') and u> a v, a'v'(q, q') are defined in Eq. I|65f) . the bare-momentum integrations corre- 
spond to the two-dimensional regions of the (fc, w)-plane where Co{q, q') — > 0, and the function g = g(fc, g') is uniquely 
defined by the parametric equation fc = lk a ^ ta > u >(q, q'). 

An example of the type of singularity l|76l) is the (5-peak which appears in the spin-singlet Coop er pair I — 1 
addition spectral function at the point (ko,u>o) — (tt, 2/i), as a result of the ^-pairing mechanism |20|. One-electron 
spectral- function singularities of the type l|77|l are studied in Ref. ^(| • It is found in that reference that some of the 
cO and si pseudofermion branch lines behave as in Eq. JZ3 for U/t -> with C c = C s = 1, where the correct C//t = 
non-interacting one-electron spectrum is recovered. The two-electron weight distribution described by expression 
(I78|l is valid when U/t — > with C c = C s = 1. Interestingly, for the dynamical structure factor it is also valid when 
U/t — > co, where such a spectral function is that of the spin- less fermions which describe the charge degrees of freedom 
in that limit. There is a qualitative difference between the spectral-function behavior described by Eq. 1(75) 1 and that 
described by Eqs. Q77f) and (|78|l : the former behavior is valid at isolated points for finite value of U/t, whereas the 
latter behavior is valid for U/t — > or U/t — > oo only. Moreover, in the case of expressions (77J| and (|78|l . as £o — * 
all the spectral weight becomes located onto the (fc, w)-plane regions where 3o(fc, oj) — > 0. 

Finally, for the n = 1 Mott-Hubbard insulator phase there are other approaches to the study of finite-energy 
spectral-weight distributions for finite values of U/t. A method for the description of the weight distribution in the 
vicinity of specific spectral features was proposed in Ref. |2l| . Also for n = 1 but for very small values of U / 1 the ID 
Hubbard model can be mapped onto the sine-Gordon model. In this limit the form factor approach was used in the 
study of the zero-momentum two-electron charge spectral- weight distribution associated with the frequency-dependent 
optical conductivity In general, the form factor approach does not apply to the ID Hubbard model. 
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VI. DISCUSSION AND CONCLUDING REMARKS 

In this paper we derived general closed-form analytical expressions for the finite-energy one- and two-electron 
spectral-weight distributions of a ID correlated metal with on-site electronic repulsion. Our results also provide 
general expressions for the one- and two-atom spectral functions of a correlated quantum system of cold fcrmionic 
atoms in a ID optical lattice with on-site atomic repulsion. In order to solve that problem, we addressed the issue of the 
dominant spectral-weight processes. That allowed the derivation of the basic weight-distribution pieces corresponding 
to surface-like, line-like, and point-like spectral features, which involved the introduction of the concepts of a border 
line and a singular and edge branch line. Our study reveals that the functional (|63|1 plays a key role in the control of 
the spectral-weight distribution. Indeed, most singularities are of power-law type and correspond to discontinuities 
of the value such a functional in the (k, a>)-plane. Moreover, that functional also controls the exceptional occurrence 
of 5-function singularities, which correspond to the (fc, w)-plane points where its value vanishes. 

For the one- and two-electron problems considered in this paper, the spectral-function expressions generated by 
processes associated with excited J-CPHS subspaces with none, one, and two active scattering centers away from the 
av = cO, si Fermi points and av ^ cO, si limiting bare-momentum values describe nearly the whole (k, w)-plane 
weight distribution. In addition to other spectral features, here we presented a detailed derivation of the spectral- 
weight distribution expressions used in the studies of Refs. 0? ]. These expressions are used in Ref. [? ] in the 
study of both the spectral-weight distributions observed for the TTF and TCNQ stacks of molecules in the quasi- 
1D organic compound TTF-TCNQ, following the preliminary predictions of Refs. 0,0,0] f° r the TCNQ spectral 
dispersions. Interestingly, the studies of Ref. [? ] find quantitative agreement with the observed spectral features for 
the whole experimental energy band width. Recently, the dynamical density matrix renormalization group method 
was used in the study of the one-electron spectral function of the ID Hubbard model j2^| for electronic densities n < 1. 
These numerical studies found spectral features consistent with those of our spectral-weight distribution expressions. 
Moreover, the results of Ref. Q are consistent with the phase diagram observed in the (TMTTF) 2 X and (TMTSF) 2 X 
series of organic compounds and explain the absence of superconductive phases in TTF-TCNQ. Other applications 
of our finite-energy spectral-weight-distribution expressions to several materials, correlated quantum systems, and 
spectral functions are in progress. This includes the use of our theoretical results in the two-atom spectral- weight 
distributions measured in ID optical lattices p|. Finally, it is shown elsewhere that the general spectral- weight 
distribution expressions introduced in this paper reproduce the low-energy behavior obtained by conformal-field theory 
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APPENDIX A: GROUND-STATE RAPIDITY FUNCTIONS 

In this Appendix we introduce simple expressions for the inverse of the ground-state rapidity functions k (q), 
A° v (q), and A°„(g) that appear in the pseudofermion energy band expressions l(T%l) - lj2U)) . These expressions are 
derived by solution of the integral equations obtained by introducing in Eqs. (13)-(16) of Ref. [j| the ground- 
state distribution functions given in Eqs. (C1)-(C3) of the same reference. The solution of these equations can be 
written in terms of the inverse functions of k°(q) and A° J/ (g) for v > 0, which we call q^ik) and q° v (A), respectively. 
Here k and A are the rapidity-momentum coordinate and rapidity coordinate, respectively. These are odd functions 
such that q®{k) = —q®.(—k) and g„„(A) = — g , ^ I/ (— A). The domains of the rapidity-momentum coordinate k and 
rapidity coordinate A are such that — tt < k < +tt and — oo < A < oo, respectively. It follows that fc°(±7r) = ±7r, 
A° (±7r) = 0, and AP av {±q^ v ) — ±oo for av ^ cO. Such relations are equivalent to q°(±ir) — ±ir, ^(±00) = ±fc^|, 
^(±00) = ±[tt — 2k f] for cv ^ cO, and g^(±oo) = ±[Ajrf — kp\\ for sv 7^ si. By suitable manipulation of Eqs. 
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(13)-(16) and (C1)-(C3) of Ref. we find, 



°+Q „ r+Q „ 

q° c (k)=k+ I dfc'$ c0 , c0 (fc',fc) ; q° su (A)= dfc'# e o,a„(A / J A) ; ^=1,2,..., (Af) 



-Q •/-Q 

g&,(A) = 2Re{arcsin^A + w?7/4fj} -J dfc' $ c0 , cv (fc', A) ; i/ = l,2,..., (A2) 

where the two-pseudofermion phase shifts $ are such that, 

- /4isinfc 4tsinfc'\ - , ¥ . - /4isinfc 4tA'\ , . . 

$ c o,co(M) = <J> c0 ,cO I — ^ — , — Jj — j; $ c0 ,^(fc,A') = $ c0 , Q „ I — J , (A3) 

£ / a ,a * ^ 4<A 4tsinfc'\ - - /4tA 4tA'\ 

$a„, c0 (A,fc') = $ Q „, c0 I — , — - — J; (A, A') = I — , J , (A4) 

and the phase shifts & a u, a'v' (p, f') ar e defined by the integral equations (A1)-(A13) of Ref. . Moreover, the 
parameters Q = k°(2kp) and B = A° 1 (fc J 7|) are such that q®(±Q) = ±2kp and q^ 1 (±B) = ±&fm and are self- 
consistently defined by the solution of the following equations, 



r+Q _ r+Q 
2k F = Q + dfc$ c0 , c0 (k,Q) ; k Fl = dk<P c0 _ sl {k,B) . (A5) 



APPENDIX B: SPECTRAL- WEIGHT DISTRIBUTIONS AND RELATED QUANTITIES 

In this Appendix we provide complementary information about the evaluation of the finite-energy spectral-weight 
distribution expressions studied in this paper. We start by expressing the generators appearing in the matrix elements 
given in Eqs. (59) and (60) of Ref. [10| in terms of the sets of canonical-momentum and bare-momentum values of 
the corresponding distribution function deviations provided in Eqs. H13|) - (|17|) . For the generators associated with 
av = cO, si pseudofermion occupancy configurations this leads to, 

F l-NF,av = LI [fl.au ff jtm ] II [e(sgn(AN^f))fl al/ + e[-sgn(ANZ F ))u, al ], (Bl) 

j = l i=l 

iFav+^qFaiy, +1 
iFau+^Fav, +1 

F J-GS,av = II F J-GS, av,t'> F J-GS, av, i = S h +1 1 II ad + -l[ 4i , J ' ( B3 ) 

F p-h,av~ F p-h,av,Li F p-h,av,i = \\_ \f\u , av fq' jo , av i fc = ±1 . (B4) 

t=±l j t =l 

For the generator corresponding to <w ^ cO, si pseudofermion occupancy configurations we find, 

4,»=fl/L; wco,si. (B5) 
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The operators I|B1|) and (|B5|I generate the elementary processes (A), whereas the operator l|B4|) generates the elemen- 
tary processes (C). The operators F_qs and F^_ GS create the whole av = cO, si pseudofermion Fermi sea, 
including the Fermi point shifts generated by the elementary processes (B), for the discrete canonical-momentum val- 
ues of the ground state and reduced-subspace excited state, respectively. In the second expression of both Eqs. I|B2(I 
and l|B3(l the products refer to the av = cO, si pseudofermions belonging to the i = sgn(q) 1 sub-branch and AqFav, l 
and Aqp al ,^ L are the bare-momentum and canonical- momentum shifts provided in Eqs. (|10|) and l|ll|) . respectively. 

The weight of the peak corresponding to the N^, av = cO, si av pseudofermion particle-hole excitation mentioned 
in Sec. IV reads, 



Aa V a = i a' av a "' ; a> a v = aa I f J '(q 1 , ■ ■ ■ , q N Ph, q i, • ■ • , q N gh) ; av = cO, si . (B6) 

Here is the lowest-peak weight given in Eq. (|48|l . is the relative weight, and qi,--- ,q N ph and 

q\, • • • , q' N ph are the corresponding set of av — cO, si pseudofermion "particle" and "hole" bare-momentum values, 

respectively. By use of the method of Ref. Q , we find that for iV^ = 1 the relative weight i|B6|l is given by, 

^•'■'-feiEJEj' °" = * sl - (B7) 

The function g(q) appearing on the right-hand side of Eq. (|B7J) reads, 



QFa 



9(13) = n 



qAa 



n 



q"i=<iF a 



ai^ = cO, si . 



(B8) 



where the correspondence between qj and is defined by the equation q~j = qj + Q®„(qj)/L with Q®„(qj)/2 given in 
Eq. |@J, qFav, ±1 and gFai^, ±1 are the Fermi points provided in Eqs. I|l()(l and (| 1 1 f> - respectively, and the corresponding 
products run over the discrete ground-state bare-momentum and excited-state canonical-momentum values, except 
for q'j — q~j and q"j — qj. Similarly, for for = 2 we find, 



, x _ 9(11)9(12) sin^(l^M) sin 2 (l2^) 



where = cO, si. The expression of the function (|B6J| is similar for increasing number of av — cO, si pseudofermion 
particle-hole processes, but since it is longer, for simplicity we do not give it here. 

We proceed by deriving the spectral-function expression l|66(l associated with creation of two active scattering 
centers away from the av — cO, si Fermi points and av ^ cO, si limiting bare-momentum values by the elementary 
processes (A). For simplicity, we consider that the two created active pseudofermion and/or hole scattering centers 
belong to the av — cO, si branches. Generalization to active scattering centers associated with av ^ cO, si branches 
is straightforward. In order to arrive to the specific form of the general i = spectral function expression l|44() for 
the excitations associated with the latter processes, we start by considering excited energy eigenstates generated by 
the elementary processes (A) whose energy AE and momentum AP relative to the initial ground state have values 
such that [u> — IAE] and [k — IAP] are small and < [oj — IAE] < and — l/v s \ < 1/v < +l/v s i. Here k and u> are 
the momentum and energy values provided in Eq. (|65|l . v is the velocity defined in Eq. I|46|l . and the corresponding 
energy AE = AE(q + Aq, q' + Aq') and momentum AP = AP(q + Aq, q' + Aq') are the general functionals given in 
Eqs. I|28l) and (|29|l . respectively, which for the present case read, 



AP = [ko + ci(q + Aq) + c[ {q 1 + Aq')} ; AP = [w + ci e av (q + Aq) + c[ e a , v .{q' + Aq')} , (BIO) 

where the bare-momentum deviations Aq and Aq' relative to the bare momentum values of Eq. (|(j5[) are small. 

The excited energy eigenstates which contribute to the weight-distribution expression i|t)6|l have momentum k and 
energy u> relative to the initial ground state. Such states are generated by the elementary processes (C) from the 
states whose momentum and energy spectra are given in Eq. (|B10ll . In order to reach the momentum k and energy 
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uj given in Eq. I|tj5fl . the elementary processes (C) generate excitations whose momentum and energy are given by 
[k — IAP] and [uj — IAE], respectively. From use of Eqs. l|H5l) and (|B10|) we find, 

[k-lAP] = -liaAq + c^Aq']; [uj - IAE] = ~l[ci v au (q) Aq + c[ v a > v >(q') Aq'} , (Bll) 
where av = cO, si and aV = cO, si. Use of these expressions leads to, 

A 1\ 7 A 771 Cl[l-?W(g')z] a / ,r , AF1 C[ [1-V av (q)z] _ 

Aq = -I [uj - IAE\ — —r ; Aq = I [uj - IAE] — — ^ ; z=l v. (B12 

For small but finite values of [uj — IAE] we can use the asymptotic expression (f61f) in the general i = spectral 
function expression Q44JI. which corresponds to the dominant contribution to the spectral function (141ft for i = 0. We 
thus find, 

B l M (k,w) « ^ 2 l c c J dq" f dq"'e(n-l[uj-lAE])e(l[uj-lAE 

F ([k- lAP]/[u - IAE]) /1[uj- lAE]yi+Co(<i,q') 
4tt^/v c0 v s1 VAtt^Vco v s i) 



l = ±l. (B13) 



Here the q" = q + Aq and q 1 " = q' + Aq' integration domains correspond to the excitation domains such that 
< [oj - IAE] < a nd -l/v sl < l/v < +l/v sl , Co = Co(?, q') is the functional JS^J, and F (z) is given in Eq. JB^J. 
By use of Eq. (|B12(I . we perform a transformation from q" and q'" to the integration variables uj' = [uj — IAE] and 
v, what leads to, 

Bj^(fc, uj) pa - - - — — / dw' / dzFo(z) ; I = ±1 . (B14) 

^ v 47r 2 C c C' s |iw(g)--tw(g , )l Jo J-i/v* W^T^IT^V 

Finally, by performing the uj' integration we reach the general expression (|66fi . 

Our next goal is the derivation of the branch-line expression defined by Eqs. I|7U|I and (|71|l . For simplicity, let us 
consider a av = cO, si branch line. The expressions correspond to k — l[ko + c± q] and uj *=s Z[oj + ci e Q j/(g)] with the 
(k, w)-plane point being located in the vicinity and just above (/ = +1) or below (I = —1) the branch line. In this 
case the initial state integrations of Eq. ( 144f t correspond to points belonging to a branch-line segment associated with 
q' = q + Aq values such that Aq is small and has both negative and positive values. The elementary processes (C) 
generate excited states of momentum k = (AP + k') and energy uj — (AE + uj') relative to the initial ground state. 
The general energy spectrum AE and momentum spectrum AP of the corresponding states of the reduced subspace 
is given in Eqs. (f28fl and (f29fl . respectively. In the present case these spectra are functions of the bare momentum 
q' = q + Aq only and read, 

AE(q + Aq) (q) + Aqa v av {q) ; AP(q + Aq) — fc + C\ q + c\ Aq . (B15) 

The velocity (f46fl is then given by, 

l(u-l[w + cxe av (q)]) 1{uj - 1[uj + ci e av (q)]) 

v = v au {q)-cx r ; Aq = -c x — , (B16) 

Aq v-v au (q) 

where we also provided the corresponding value of Aq as a function of v. 

For sgn(w) = — sgn(g) one has that v e (— oo, — v s i) if sgn(g) 1 = +1 and v € (v s i, oo) when sgn(g) 1 = — 1. 
By use of Eq. (fB16f) we then find that the limiting values — sgn(g) oo and — sgn(q) v„i correspond to Aq = and 
Aq = sgn(g) c\1(uj — 1[ujq + c\ e au (q)])/[v s i + \v au (q)\], respectively. In turn, for sgn(t>) = sgn(g) there are two 
cases. When |tw(?)| < v si [1 — l(w — ^[ w + c i e av(l)])/^\ we nn d that v € (—oo, —v s i) for sgn(g) 1 = — 1 and 
v € (v 3 i, oo) for sgn(g) 1 = 1. In contrast, if |tw(<7)| > v 8 i [1 — l(u> — 1[ujq + ci e av (q)] )/0] we find that v £ 
(-co, -[v av (q)[/[l-l(uj - l[uj Q + cx e av (q)])/ft]) for sgn(g) 1 = -1 and v £ (\v av (q)[/[l-l(uj-l[uJo + ci e av (q)])/CL], oo) 
for sgn(g) 1 = 1. The limiting values sgn(g) oo, sgn(g) v 8 i, and v av (q)/[l — l(u> — 1[ujq + ci e al/ (q)])/fi] correspond to 
Aq = 0, Aq = -sgn(g) c t l(u) - 1[uj + c t e au (q)])/[v sl - \v av (q)[], and Aq = -c t [Q - l(u> - 1[uj + c t e au (q)])]/v au (q), 
respectively. By use of the general spectral function expressions (I4T) and (|4^|) for i = we then arrive to, 
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sgn(g) ci 
2ir C r C, 



M 



v s \ [1 - J - \v av (q)\ 



<j+sgn((j) c_ 



»«l + l«oti/(9)IJ 



q— sgn(g) ci — L - 1 ^ — i <£&£*±L> 



dq' 



v,l-lvcv(.q)l) 



+ 0(\v av (q)\-v s i [1 - 



l(u; - 1[ujo + ci e my (<?)]) 
O 



+sgn(<j) c: 



x B' w - Z[w + ci e aI/ (g) + Aqci «„„(?)], tw(<?) - ci 



/[cj - l[u + ci e at ,(g)] 
A 9 



; l = ±l. (B17) 



Next, we change the integration variable from q' to v and reach the following spectral-function expression, 
, sgn(g)l r f-^)^ 

/7TO c O s k J-sgn( g ) oo 



+ 0(u a i[l j-K^wl 



'( - V s i [1 - 



?(w - 1[uj + a e av (q)}) 



9. 



sgn(g) oo 

sgn(g) w,i 
sgn(g) oo 



dv 



dv 



1 



x B 1 ' (uj - 1[uj q + Cl e a „(q)}) 



v - v av {q) 



l = ±l. 



(v - v av (q)) 2 



(Bit 



Changing the integration variable from v to z — 1/v and using of the asymptotic expression 1(611 of Eq. (|B18|) for 
small but finite values of (u — 1[uiq + c\ e av (q)]), leads finally to expressions 1(70(1 and (|71|l . (The theta functions added 
to the expression (|70(l define the energy range where it is valid.) Similar procedures lead to these expressions for the 
av ^ cO, si branches provided one uses the corresponding parametric equations given in Eq. 1(69(1 . 

The derivation of expressions ((73(1 and l|74(l . which refer to k and u> values in the vicinity of the end points of a 
av = cO, si branch line such that w « iv av (k — Iko) + lu>o, follows the same steps as for the general branch-line 
spectral-function expression. When Aw is small but finite the convolution function (|72|l can be written as follows, 



av = cO, si; t = ±1; / = ±1 , (B19) 



where we used the relation [ZA<fti ne /27r] 2A <*" = 1/N a . Here the functional £° L is given by, 

CL, t = 2A^ + 2A+i + 2Ari = C°„ - 2A^ ; av = cO, si ; u = ±1 , 
and the function t (z) reads, 



(B20) 



42A - ; [ n n r(2A i ', .) 

1 -1 i'=±l v o!v<) 



ct'v' — cO, s 



ix { n 



-i if! ii v av 
1 — X + I V av Z — U X 



dx 

V'=±l 
2Al'-,-l 



fc) ( 1 — X + i v av z — Li X 

Vrvv 



2x 



2A1„-1 



(B21) 



Furthermore, use of the same procedures as above but with the use of the functions (|72|l . ((B19(l . and (|B2 1|) instead of 
(I45|l . 161(1 . and ((62(1 . respectively, leads straightforwardly to expressions 1(73(1 and ((74(1 . The main difference is that the 
factor l/Na which multiplies the q' integration factor L/2n to give 1/2tt on the right-hand side of Eq. ((B17(l arises 
from the state summation of Eq. 144|) , whereas in the case of the constant 11741 it comes from expression ((B191 . 
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